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ABSTRACT 

The far-infrared (FIR) and radio luminosities of star-forming galaxies are linearly correlated over a 
very wide range in star formation rate, from normal spirals like the Milky Way to the most intense 
starbursts. Using one-zone models of cosmic ray (CR) injection, cooling, and escape in star-forming 
galaxies, we attempt to reproduce the observed FIR-radio correlation (FRC) over its entire span. The 
normalization and linearity of the FRC, together with constraints on the CR population in the Milky 
Way, have strong implications for the CR and magnetic energy densities in star-forming galaxies. We 
show that for consistency with the FRC, ~2% of the kinetic energy from supernova explosions must 
go into high energy primary CR electrons and that ^10% - 20% must go into high energy primary CR 
protons. Secondary electrons and positrons are likely comparable to or dominate primary electrons in 
dense starburst galaxies. We discuss the implications of our models for the magnetic field strengths 
of starbursts, the detectability of starbursts by Fermi, and CR feedback. Overall, our models indicate 
that both CR protons and electrons escape from low surface density galaxies, but lose most of their 
energy before escaping dense starbursts. The FRC is caused by a combination of the efficient cooling of 
CR electrons (calorimetry) in starbursts and a conspiracy of several factors. For lower surface density 
galaxies, the decreasing radio emission caused by CR escape is balanced by the decreasing FIR emission 
caused by the low effective UV dust opacity. In starbursts, bremsstrahlung, ionization, and Inverse 
Compton cooling decrease the radio emission, but they are countered by secondary electrons/positrons 
and the dependence of synchrotron frequency on energy, which both increase the radio emission. Our 
conclusions hold for a broad range of variations on our fiducial model, such as those including winds, 
different magnetic field strengths, and different diffusive escape times. 

Subject headings: cosmic rays - infrared: galaxies - galaxies: magnetic fields - galaxies: starburst - 
gamma rays: galaxies - gamma rays: general - radio continuum: galaxies 



1. INTRODUCTION 

The far-infrared (FIR) and radio luminosities of star- 
forming galaxies lie on a tig ht empirical relatio n, the 



FIR-radio co r relati o n" (FRC: Ivan der Kruitlll97li 
de Jong et ail 119851: iHelou et al.l 119851: ICondonl 



1973 



1992 



Yun et al.1 12001| ). The FRC spans over three decades 



in luminosity, remaining roughly linear across the range 
lO 9 ^© < L < 10 12 5 Lq, from dwarf galaxies to local 
ultra-luminous in frared galaxies (ULIRGs) like Arp 220 
(|Yun et al.ll2001| ). At low luminosities (L < 1 9 L^), the 
corre l ation shows evidence of non-l inearity (|Yun et al.l 
120011: [Be!i|[2003t iBeswick et al.ll2008t ). The galaxies that 
make up the FRC span a large dynamic range, not just 
in bolometric luminosity, but also in gas surface den- 
sity 4 (0.001 g cm~ 2 < S g < 10 g cm~ 2 ), photon energy 
density, and presumably magnetic field strengt h. From 
the observed Schmidt law of star formation ([Schmidt] 
[1951: iKennicuttl [1998), the range in gas surface den- 
sity corresponds to a range of at least 4 x 10 5 in pho- 
ton energy density. Not only does the FRC hold on 
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galactic scales, but it exists for regions within star- 
forming galaxies down to a few hundred parsecs (e.g.. 
Beck fc Gollal [19881: iBicav fc Heloul [l990l: i Murphv et al.l 
2006al iPaladino et al.ll2006l : iMurphv et al.ll2006b[ T20081 . 



Star formation drives the FRC. Young massive stars 
produce ultraviolet (UV) light, which is easily absorbed 
by dust grains. The dust reradiates in the FIR, produc- 
ing a linear correlation between star formation rate and 
the FIR luminosity, if the dust is optically thick to the 
UV light. The non-thermal GHz radio continuum emis- 
sion observed from star-forming galaxies is synchrotron 
radiation from cosmic ray (CR) electrons and positrons, 
believed to be accelerated in supernova (SN) remnants. 
Since SNe mainly occur in young stellar populations, this 
means that star formation is directly linked to normal 
(non-active ga lactic nucleus) radio emission (reviewed in 
lCondod[l992h . 

In this paper, we model the FRC, over its range in 
physical parameters from normal star-forming galaxies 
to the densest and most luminous starbursts. Our moti- 
vation is that the normalization and linearity of the FRC 
has strong implications for the physical properties of star- 
forming galaxies and the CRs they contain. For example, 
we can use the radio emission to estimate the energy in- 
jection rate and equilibrium energy density of both CR 
electrons and protons. This is important because the 
CR pressure is kno wn to be dynamically important in 
the Milky Way fe.g iBoulares fc Coxlll990D . and possibly 
starburst galaxies (S ocrates et al.l[2008T l. Furthermore, 
we can use the inferred CR proton energy density to cal- 
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culate the flux of gamma-rays from pion produc tion in 
the galaxies' host interstellar medium (ISM; e .g., iTorresI 
12004 iThompson. Quataert. fe Waxmanll2007h . Finally, 
the radio emission also constr ains the magnetic fiel d 
strength in galaxies on the FRC (j Thompson eTaD[2006). 

Finding the causes of the linearity and span of the 
FRC is the other main purpose of this paper. The FRC 
is affected by the density of CRs and the environment 
they propagate through. For the Milky Way, the prop- 
agation of CRs has been w ell studied, both observation - 
ally and theoretically (e.g.. IStrong fc Moskalenkolll998l ). 
However, given the vast range of environments in star- 
forming galaxies, it is not clear that our knowledge of CR 
propagation in the Galaxy can be extrapolated across 
the entire FRC. Therefore, one aspect of our task in ex- 
plaining the FRC is determining the extent to which the 
properties of CR injection, such as the initial spectral 
slope and proton-to-electron ratio, and CR propagation, 
such as the rate of escape by diffusion, can apply to all 
star-forming galaxies. 

The diversity of star-forming galaxies on the FRC and 
the tightness of the correlation may imply a deeper, sim- 
pler princi ple at work. In the calorimeter theory first 
proposed in lVolkl (|1989l) , the CR electrons lose all of their 
energy before escaping galaxies, with most of the energy 
radiated as synchrotron radio emission. Thus, galaxies 
are electron calorimeters, with the energy in CR elec- 
trons being converted into an observable form. Calorimc- 
try also requires that galaxies on the FRC are optically 
thick to UV light from young stars, which is reradiated 
in the FIR. These galaxies would therefore also have to 
be UV calorimeters. If both electron calorimetry and UV 
calorimetry hold, and if synchrotron is the main energy 
loss mechanism, then the ratio of FIR to radio emission 
is simply the ratio of total starlight produced to the to- 
tal energy supplied to CR electrons, which is naively ex- 
pected to be a constant fraction of the energy from SNe, 
accounting for the FRC. 

Calorimeter theory has been questioned, however, both 
in its assumptions and its implications. For example, 
the assumption that all normal galaxies are optically 
thick to UV light is probably false: the observed UV 
luminosity of normal star-forming galaxies is compara- 
ble to the obse rved FIR lumino s ity at low overall lu- 
minos ities (e.g-IXu fc Buatlll995t lBeill2q03t iBuat et al.l 
I2005t iMartin et al.l 120051; iPopescu et al.l 12005ft . Nor is 
electron calorimetry believed to hold in the Milky Way 
(and presumably similar galaxies) , since the inferred dif- 
fusive escape time is shorter than the typical estimated 
synchrotron coolin g time (see equat ions [5] and [T2l later in 
this paper; or, e.g.. [Lisenfeld et aTlll996af) . 

Even in cases when calorimetry holds, the implica- 
tions of standard calorimeter theory may conflict with 
observations. A long-standing problem with the pre- 
dictions of calorimetry has been the radio spectral in- 
dices of star-forming galaxies. If electron calorimetry 
holds, then the synchrotron cooling timescale is much 
less than the escape timescale. The electron popula- 
tion will then be strongly cooled with a steep spec- 
trum. For an initial injection spectrum of N cx E~ p 
where p « 2 — 2.5 and a final synchrotron-cooled steady- 
state spectrum N oc E~ v , this would imply a syn- 
chrotron spectrum of F v oc v~ a with a spectral index 
of a = (V - l)/2 = p/2 w 1.0 - 1.2. The observed spec- 



tral indices are 0.7 — 0.8 for normal galaxies, suggest- 
ing that, contrary to calori meter theory, ele c trons e scape 
before losing their energy. iLisenfeld et al.l (|1996a|) con- 
sider a modified calorimeter model for normal galaxies 
that includes escape co mparable to cooling losses, and 
ILisenfeld fc Volkl (|2000[ ) suggest that SN remnants in the 
galaxies can flatten the observed radio spectrum. 

More drastically, s everal non-calorimeter theories have 
been proposed (e.g.. lHelou fc Bicavl fl993: Nik las fc Beckl 
1997), often involving a "conspiracy" to maintain the 
tightness of the FRC. A potential pitfall of non- 
calorimeter models stems from the enormous dynamic 
range in physical properties for galaxies on the FRC. 
For example, inverse Compton cooling alone is very 
quick in starbursts, implying that electrons cannot es- 
cape from these galaxies before losing most of t heir en- 
ergy (jCondon et al.lll991l ; IThompson et al.ll2006D . 

Typical explanations of the FRC leave out two under- 
appreciated but important effects: proton losses and non- 
synchrotron cooling. Models of individual starbursts, 
which have gas densities 10 3 — 10 4 times higher than 
the Milky Way, predict that CR protons lose most of 
their energy to pion creation as they interact with the 
ISM. When a CR proton collides with a proton in the 
ISM, it produces a pion, either charged (ir + or or 
uncharged (tt°). Neutral pions decay into gamma rays, 
so that pion losses should act as a source of gamma-ray 
luminosity in starbursts. Charged pions ultimately de- 
cay into neutrinos (which may eventually be observed 
with neutrino telescopes), as well as secondary electrons 
and positrons. Therefore, dense starburst galaxies are 
expected to be proton calorimeters: essentially all the in- 
jected energy in CR protons ends up converted to gamma 
rays, neutrinos, and secondary electrons and positrons. 
Proton calorimetry would also imply that, unlike the 
Milky Way, secondary electrons and positrons may dom- 
inate over primary electrons and positrons , depending 
on th e ratio of injected protons to electrons ([Rengaraianl 
120051 ; secondary electrons and positrons are found to be 
more abundant than primary electrons in starburs ts by 
Torres! 12004 iDomingp-Santamaria fc Torresl I2005L and 
de Cea del Pozo et al.l 120091 ) . Because secondary elec- 
trons and positrons radiate synchrotron, their presence 
poses a problem for any explanation of the FRC that 
requires the CR electron density to be directly propor- 
tional to the star formation rate, includ ing both stan- 
dard c alorimeter theory and the theory of Ni klas fc Beckl 
(fl997h . 

On the other hand, bremsstrahlung, ionization, and 
IC may all be mor e imp ortant in starburst galax- 
ies. IThompson et al.l (|2006f ) point out that cooling by 
bremsstrahlung and ionization tends to flatten the radio 
spectra, thus saving calorimeter theory from the spectral 
index argument, at least for starbursts. However, the 
energy CR electrons lose to bremsstrahlung, ionization, 
and IC cannot go into synchrotron radio emission, an ob- 
stacle for any theory that assumes radio emission is di- 
rectly proportional to the injected power of primary CR 
electrons. Therefore, even electron calorimetry and UV 
calorimetry are not enough to guarantee a linear FRC. 

We address these issues with one-zone numerical mod- 
els of CRs in star-forming galaxies. These models include 
CR escape as well as the main cooling processes and sec- 
ondary production, a combination that has not to our 
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knowledge been done over the entire span of the FRC. 
CRs in individual galaxies have been studied with similar 
one-zone models that fit the emissio n across the electro- 
magnetic spectrum (e.g., Ar p 220 in lTorresll2004t M82 in 
Ide Cea del Pozo et al.1 12009) . but in this paper our focus 
is on the FRC itself and not any individual galaxy. Our 
one-zone approach allows us to efficiently parameterize 
unknown quantities like the magnetic field strength, and 
to try a large number of scenarios. The primary inde- 
pendent variable in our calculations is the gas surface 
density S ff , which controls both the photon energy den- 
sity through the observed Schmidt law and the average 
gas density. These simplifying parameterizations allow 
us to qualitatively understand the FRC over the range 
of star- forming galaxies, although it ignores deviations 
and complications that may be important for individual 
galaxies. 

We first describe the calculations necessary to find 
the CR spectra and observables for each galaxy (Sec- 
tion [2|). We review the effects of each parameter on 
the observables (Section [3]), before presenting our results 
(Section [4j. We discuss the implications of our work 
for CR physics (Section [5]), including whether calorime- 
try is correct (Section l5.ip . what causes the FRC (Sec- 
tion 15.21) , predictions for the FRC at other frequencies 
(Section 15 . 3|) . the spectral slope problem (Section I5.4[) . 
the gamma-ray luminosities of galaxies (Section l5.5[) . and 
whether CR pressure and magnetic pressure are impor- 
tant as feedback mechanisms in galaxies (Section I5.6j) . 
We finally summarize our results (Section |6]). In Ap- 
pendix we present results from a suite of variants on 
our standard model, and show that those consistent with 
the FRC have similar parameters to our standard model. 
For the reader's convenience, we list the symbols we use 
in our calculations and discussions in Table [TJ 

2. PROCEDURE 

We construct one-zone leaky box models of galaxies 
across the dynamic range of the observed FRC. We treat 
star-forming galaxies as homogeneous disks of gas, char- 
acterized by a column density E 9 , a star formation rate 
surface density Ssfr, and a scale height h. We solve 
the steady-state diffusion-loss equation for the equilib- 
rium CR spectra of primary and secondary electrons and 
positrons, as well as primary CR protons. 5 

Under these simplifying assumptions, the diffusion-loss 
equation for CRs becomes 

where E is the total energy, N(E) is the CR spectrum, 
tlife(E) is the energy-dependent lifetime to diffusive or 

5 In the steady-state approximation, the dN/dt term in the 
diffusion-loss equation is assumed to be small. For the Milky Way 
(and presumably other normal spirals), the CR flux is known to be 
constant within a factor of ~ 2 for the last billion years from solar 
system studies (Arnold ct al. 1961; Schacffcr 1975). The tightness 
of the FRC combined with the long timescales for galactic evolution 
also imply the CR population in normal galaxies is steady state. 
Additionally, in extreme star bursts, the IC cooling time alone for 
GHz-emitting CR electrons (< 10 4 yr) is much shorter than the 
characteristic timescale for the system to evolve. Therefore, we ex- 
pect that the steady-state assumption is valid. However, we note 
that in weaker starbursts, where the cooling and escape times for 
CRs are several Myr, evolution m ay be important and t he steady- 
state approximation may fail (sec Liscnfcld et al. 1996b). 



advective escape from the system, Q(E) is the CR source 
term, and b(E) = —(dE/dt) is the rate of energy loss for 
each particle. The equilibrium CR spectrum is a compe- 
tition at every energy between injection, cooling, and es- 
cape losses. If the injected CRs initially have a spectrum 
of the form Q(E) oc E~ v and if escape is insignificant 
(tm c (E) — > oo), the final spectrum will have the form 
N(E) oc E l ~ v /b(E). If instead cooling is insignificant 
compared to escape (b(E) —> 0), the final spectrum will 
have the form N(E) cx E~Pt mc {E). 

We solve the general form of equation ([T]) numeri- 
cally using a Green's fun ction for CR protons, elec- 
trons, and positrons (see iTorred [2004). We include 
sync hrotron, IC, bremsstrahlung, and ionizat ion losses 
(e.g.. iRvbicki fc Lightmanlll979t lLongairlll994D . For CR 
protons, we also include pion losses in t\a c {E) due to 
in elastic p roton -proton collisions using the formalism 
of iTorresI (l200l! 6 The pub l icly available GALPROP 
code 7 (tStrone :k Moska lenkcj 119981 : IStrong et all 120001 : 
Moskalcn ko et al.1 12002) is used to calculate the differ- 
ential cross section for electron or positron production 
from proton-proton collisions, as well as for calculating 
the spectrum of 7-rays produced by the decay of sec- 
ondary 7T° mesons. We also include knock-off electrons 
from CR p roton collisions with atoms in the ISM (see 
lTorresll2004h . 

2.1. Primary CR Injection Rates 

We assume that both the primary CR electrons and 
protons are injected into galaxies with power law spectra 
Q{E) = CE-v with 1 < 7 < 10 6 (where 7 = E/{mc 2 ) 
is the Lorentz factor), and we consider initial spectral 
slopes p in the range 2.0 < p < 2.6. Integrating the in- 
jection spectrum times the kinetic energy per particle K 
gives the total power injected per unit volume for each 
primary species, ecR — J KQ(E)dE, to set the normal- 
ization. Energetic and escape losses produce the final, 
steady-state spectrum as determined by the solution to 
equation (fT]). 

In order to normalize the CR injection spectra, we as- 
sume that a constant fraction £ and n of the kinetic en- 
ergy of SN explosions (E 5 i — i?sN/10 51 erg) goes into 
accelerating primary CR electrons and protons, respec- 
tively. The CR electron and proton emissivities ccr are 
then proportional to the emissivity in starlight photons, 
e p h, produced by star formation, when averaged over 
the star formation episode. 8 Following th e discussion in 
iThompson. Quataert. fc Waxmanl ((2007) , we calculate 
the starlight emissivity (here, in units of erg s _1 cm~ 3 ) 
as 

e ph = £E SFR c 2 (2/i)- 1 , (2) 

6 Including pion losses in the cooling term b(E) instead is for- 
mally incorrect, because the losses are catastrophic instead of con- 
tinuous. For p = 2.0, the resulting proton spectra are nearly iden- 
tical, but for p = 2.6, including pion losses in b(E) decreases the 
proton spectrum N(E) by ~ 40% when pion losses are strong. 

7 Specifically, we use the "PP_MESON" subroutine, which 
calculates the cross sections for electron and photon produc- 
tion through pion p roduction. GALPROP is available at 
http:/ /galprop. stanford.edu 

Though we assume that the SN rate is proportional to the 
starlight, in reality, there will be a lag between the first massive 
stars and the first SNe when there will be very few CRs (e.g., 
IRoussel et"aL|[200l) . which we do not account for. 
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where e = 3.8 x 10~ 4 is a dimensionless initial mass func- 
tion (IMF) dependent constant that relates the luminos- 
ity in young stars to the instantaneous star formation 
rate, c is the speed of light, h is the CR scale height (in 
cm), and Ssfr is the surfa ce density of star formation in 
cgs units of g s cm" " 2 (|Kennicutdll998| ). Equation [2] 
essentially says that some proportion of the mass that 
forms stars is converted into starlight. We take the sur- 
face density of star formation Esfr from the observed 
Schmidt law ()Kennicuttlll998D . 
The cmissivity in CR electrons can then be written as 



e C R, = 9.2 x 10 5 E 51 i[j ir e ph 



0.01 



(3) 



where Vi7 = (r SN /e)/(17M~ 1 ) and r SN is the SN rate 
per unit star formation. Similarly, the total emissivity of 
the primary CR protons can be written as 

ccr, p = S ecR, e 

= 9-2xlO-^ 17 e ph (A)(X), (4 ) 

where 6 = rj/£ is the ratio of the total energy injected 
in CR protons to that in CR electrons per SN. Although 
we have normalized £ = 0.01 and 5 = 10 in the above 
expressions for reference, one purpose of this paper is 
to show explicitly that numbers in this range are in fact 
compatible with observations of radio emission from star- 
forming galaxies. 

2.2. Environmental Conditions 
2.2.1. Escape 

The CR lifetime (tufe) for both primary electrons and 
protons in equation ([1} is uncertain, and probably varies 
from normal sp irals like our ow n, where losses are mainly 
diffusive (e.g., lLongairl |1994[ ). to dense starbursts like 
M82 and Arp 220, where CRs ar e likely advected in 
a lar ge-scale galactic wind (e.g., iSeaauist fc Odegardl 
[1991. 

We use a prescription for diffusive losses motivated by 
observations of beryllium isotope ratios ("CR clocks") 
at the Solar Circle, which suggest that the confine- 
ment timescale for CR protons with E > 3 GeV is 
(e.g.. | Garcia-Munoz et aT] Il977t lEngelmann et all 119901 : 
IConnelllll998l:lWebber et al.ll2003ft 



t dm {E) = 26 Myr 



E 



3 GeV 



-1/2 



(5) 



Identifying t\if e with a diffusion timescale on a physical 
scale of ~ kpc implies that CR protons diffuse in the 
ISM of the Galaxy with a scattering mean free path of 
order ~ pc. Although the behavior of t\i( e at lower en- 
ergies is uncertain because of the effects of solar modu- 
lation (compare, e.g., Engelmann et al. 1990 & Webber 
et al. 2003), we use equation §5§ for all fiducial mod- 
els employing diffusive losses. Variations to the diffusion 
constant are considered in Appendices IA.6IIA.8l 

There is ample evidence for large-scale mass-loaded 
winds in starb urst g alaxies wi th £ > 0.05 g cm~ 2 
(|Heckman et al.ll2000UHeckmanll2003l ). These winds can 



advect CRs out of their host galaxies on a short timescale 
with respect to equation (0 , thus affecting both the pre- 
dicted emission and the overall equilibrium energy den- 
sity of CRs. For this reason, we consider in Section IA.2I 
models of starburst galaxies with 



t 



wind 



^wind 



300 kyr (E fl > 0.05 g cnT 



(6) 



Note that 1 g s 



1.5 X 10 11 M Q yr" 1 kpcT 



where we have taken h = 100 pc, a wind speed of 
v w ind = 300 km s _1 , and a cutoff between starburst and 
non-starburst galaxies of S ff > 0.05 g cm -2 as reference 
values. The combined escape time in equation (JT]) is then 
given by = + t^ nd . 

2.2.2. Scale Height 

In our one-zone models, the galaxy scale height repre- 
sents the volume in which the CRs are confined for tuf e 
(eq. [TJ . Because we specify the properties of galaxies by 
their surface density, h is also important in determining 
the average gas density seen by CRs, which, in turn, is 
important for bremsstrahlung, ionization, and pion losses 
(see Section [2X3]). 

We adopt h — 1 kpc for normal galaxies (S s < 
0.05 g cm -2 ), altho ugh we consider several other scale 
heights in Appendix IA.3I There are several relevant scale 
heights which are not identical, and we must choose one 
for our one-zone model. CRs are injected in the gas disk 
(h fa 100 pc), but we do not use this scale height since the 
CRs diffuse and emit synchrotron outside the gas disk. 
The observed beryllium isotope ratios, as interpreted by 
CR diffusion models imply that the C Rs have a scale 
height of 2 - 5 kpc in normal galaxies (ILukasiak et alJ 
[1991 [Webber fc Soutouj[l99l [Strong et all 120001) . The 
magnetic fi eld scale heights of normal g alaxies are also 
several kpc (|Han fc Qiao 1994; Beck 2009). Finally, radio 
emission in most normal galaxies come from two disks: a 
thin disk wit h h ~ 0.3 kpc and a thick disk with h sa 
2 kp c (e.g., [B eucrmai m et alJ 119851 : iDumke fc Krausel 
119981 : iHeesen et aU 120091 ). On average, the thin and 
thick disks emit the same radio power at ~ 1.4 GHz, 
and one component fits of the radio emission of normal 
galaxies general ly find a one-component scale height of 
~ 1 - 1.5 kpc ([Dumke et al.lll995l l2000t iKrause et afl 
[20061) . 

Since we are most interested in the radio emission of 
normal galaxies to explain the FRC, we use the value 
of h = 1 kpc from the one component fits. However, 
this one-zone approach does not capture all the relevant 
physics of the CRs. In particular, the escape time in 
Equation © applies to the entire CR halo. Escape from 
the radio-emitting regions is likely quicker and is prob- 
ably underestimated in our models. Conversely, a vari- 
ant on our fiducial model with large h in Appendix I A. 3 1 
probably overestimates the synchrotron losses in nor- 
mal galaxies, since it does not account for the lower 
magnetic field strengths far from the midplane. Note 
that h = 1 kpc implies a vertical diffusion constant of 
D z w 7 x 10 27 cm 2 s^iE/GcV) 1 / 2 (for 1.4 GHz emitting 
electrons in the Milky Way, D z 1.3 x 10 28 cm 2 s" 1 ; 
compare with the valu es in iDahlem et al.l 119951 and 
iPtuskin fc Soutoull[r998h . 

Starburst galaxies (S ff > 0.05 g cm -2 ) are considerably 
more compact, both in terms of their star forming regions 
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and in terms of their CR c onfinement zone, and for t hem 
we adopt h = 100 pc (e.g.. iDownes fc Solomonl[T998h . 

2.2.3. ISM Density 

Estimates from beryllium isotopes imply that the av- 
erage density experienced by CRs is about one fifth to 
one tenth that of the Galactic disk. This may be because 
the ISM is clumpy and the CRs avoid the clumps, or be- 
cause the CRs spend significant time in the low-density 
Galactic halo. 

CRs do not necessarily travel through gas with the 
mean ISM density. The actual average density CRs ex- 
perience depends upon the injection and propagation of 
the CRs, which depends on the small-scale ISM struc- 
ture in galaxies and starbursts. For example, most of 
the volume of the ISM in galaxies is low density mate- 
rial, and we can imagine the CRs are injected into this 
low density phase and lose their energy before encounter- 
ing high density clumps. Then the density experienced 
by CRs is lower than the average gas density. Conversely, 
we can imagine that CRs are preferentially injected into 
high density clumps, and are confined there by magnetic 
fields in the clumps, in which case, the CRs experience a 
higher density than the average density of the galaxy or 
starburst. 

For this reason we include a parameter /, to account 
for these unknown propagation and injection effects, de- 
fined by 

n cS = f{n), (7) 

that measures the effective density "seen" by CRs (n e g) 
with respect to the average density of the CR confine- 
ment volume, (n) = T, g /(2h). For / > 1 or / < 1, the 
CRs traverse over- or under-dense material compared to 
(n), respectively. Note that we are defining / with re- 
spect to the CR confinement volume and not the gas 
disk. 

The primary importance of the parameter / is in deter- 
mining the importance of bremsstrahlung and ionization 
losses for CR electrons and positrons, and of pion pro- 
duction from inelastic proton-proton collisions. 

Note that even though both the star formation rate 
and the magnetic fields of galaxies are taken to depend 
on the surface density (see Section 12.2. 5|) , they are as- 
sumed to be independent of /. This means that the 
radiation energy density and the magnetic field the CRs 
experience are assumed to be average, while the CRs are 
allowed to traverse through underdense or overdense ma- 
terial. Although we allow ourselves this freedom in the 
modeling, it turns out that the models with / ~ 1 are 
most consistent with observations in the Milky Way. For 
example, our adopted gas surfa ce density at the Sola r 
Circle, E g = 2.5 x 10~ 3 g cm- 2 (|Boulares fc Coxlll990D . 
and scale height h = 1 kpc imply an average number den- 
sity of (n) — 0.24 cm -3 . Since the CRs are inferred to 
trave l through material of density n e s ~ 0.2 — 0.5cm~ 3 
(e.g.. IConnelllll998t ISdilickeiserll2002[ ). this implies that 
jfwl. 

2.2.4. Interstellar Radiation Field 

The interstellar radiation field is important for deter- 
mining the IC losses for CR electrons and positrons. 
The primary contributions to the interstellar radiation 
field are starlight and the cosmic microwave background 



(CMB). The latter is particularly important for low sur- 
face brightness galaxies where it dominates starlight. 
Both sources of radiation are included in all models. 

When the galaxy is optically-thin to the re-radiated 
FIR emission from young stars, then the energy density 
in starlight, which dictates the IC cooling timescale, is 
simply 

£W = i^/c = £ S SFR c (8) 

/ E \ 1A 

= 3xl0~ 9 ^ erg cm" 3 , (9) 

\gcm~ 2 J 

where the surface density of star formation Eg FR is con- 
nected to the average gas surface density by the Schmidt 
law. For large gas surface densities (E s > 0.1 — 1 g cm~ 2 ) 
galaxies become optically thick to the reradiated FIR 
emission and 

C/ph,* = (tfir + l)F*/c = (r Fm + l)eE SFR c (10) 

/ E \ 1A 

= 3 x 10- 9 (r Fm + 1) a —, erg cm^ll) 

\g cm z / 

where t F ir = k F irE 9 /2 is the vertical optical depth, and 
K F rR is the Rosseland mean dust opacity. For parame- 
ters typical of starbursts and ULIRGs, k F ir 1 — 10 cm 2 
g" 1 for Galactic dust -to-gas ratio and solar metallicity 
(|Semenov et al. I [20031 ). For our standard models (Sec- 
tion SH]), we assume that the CRs are always in opti- 
cally thin regions, so that equation ([9]) hol ds. H owever, 
we discuss models with t F ir > in Section IA.4I 

2.2.5. Magnetic Fields 

A primary motivation for this work is to determine how 
the average magnetic energy density of galaxies scales 
from normal galaxies like our own to dense ULIRGs 
like Arp 220. Observations of Zeeman splitting in 
ULIRGs supports a relatively strong sc aling of magnetic 
field strength with gas surface density (fRobinsha w et all 
|2008| ) . To test a suite of models for consistency with ob- 
servations, we parametrize the global average magnetic 
field of galaxies as 

^ g - 6 ( o.0025 S g 9 cm- 2 ) Q ^ < 12 > 
where a is determined from comparing with the FRC, 
and where the normalization has been chosen to 
match fiducial numb e rs at the Solar Circle (as i n 
IBoulares fc Coxl H^90t IStrong et all l2000l lBecll200l . 
The magnetic field energy density is then just Ub = 
B 2 /(8tt). The E^ dependence is motivated by the Parker 
instability: the magnetic energy density cannot exceed 
the gas disk midplane pressure 7rGE 2 , or else the mag- 
netic field will buoy up out of the disk and escape (|Parkerl 
1966). A natural scaling for B given the Parker limit 
would be B oc E fl . The E° scaling also arises if the mag- 
netic field is in equipartition with the starlight, because 
the Schmidt law implies that U p h oc E^ 4 ; in this case 
a = 0.7. We consider 0.4 < a < 1.0. We assume that the 
magnetic field is constant within the confinement volume 
of scale height h, a reasonable assumption based on the 
observed radio halos of galaxies an d Galactic pulsar rota- 
tion measures (|Han fc Qiaolll994l ). We also consider two 
other parameterizations of the magnetic field, B oc p a 
and Ub = U p h, in Section |A~T1 and Section [A~5l respec- 
tively. 
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2.3. Observables and Constraints 

We use two broad conditions to select successful mod- 
els. We first ask if the model satisfies the FRC. However, 
since the FRC alone does not necessarily demand CR 
protons at all, a second constraint is needed to fix the 
overall CR proton normalization. We consider two sets 
of constraints on the protons, either using Earth-based 
measurements of CRs, or observations of the entire Milky 
Way. 

1. Reproduce the FIR-radio correlation. The nonther- 
mal radio luminosity is calculated directly from our 
synchrotron spectrum as ve v at v = 1.4 GHz. We 
do not include the thermal free-free contribution 
to the radio luminosity; however, the thermal ra- 
dio luminosity of most galaxies is typically small 
at GHz frequencies. Nor do we consider the effects 
of free-free absorption. The total infrared (TIR) 
luminosity 10 is e p h[l — (1 — exp(— ttjv))A"uv]> with 
the UV optical depth t\jv through the entire disk 
calculated as KuvEg- We adopt a UV opacity of 
500 cm 2 g , which is roughly appropriate at wave- 
lengths of - 1000 A (T off « 30, 000 K ] and Galac- 
tic m e tallicity an d dust-to-gas ratios ([Li fe Draind 
1200 It iBel 120031 use k uv = 190 g cm" 2 , using 
a smaller dust-to-gas ratio and 1500 A). Then, 
£'rm/-£'radio is simply the ratio of these luminosi- 
ties, and can easily be converted into qvm, an ob- 
servable quantity we calculate 11 as 

/ £tir \ 



9fir = log 



10 



- I - 3.67 (13) 

\ -^radio / 

and defined in lHelou et al.l (|1985h . The normaliza- 
tion o f the FRC is L T m/L Tll dio = 9x 10 5 (|Yun et all 
2001), which we match by adjusting £ appropri- 
ately, therefore fixing the primary CR electron in- 
jection rate in galaxies (Section |2~T]) . 

Once we have the ratio LTiR/^radio, our primary 
constraint is that we require a linear FRC to exist. 

We require that 

max(L T IR/£radio) 



min(L T i R /L radio ) 



< 2. 



(14) 



2. Fix the proton normalization. We then use two 
sets of constraints to fix the proton normaliza- 
tion, the local constraints and the integrated con- 
straints. Each is considered independently for each 
model. For simplicity, the proton normalization is 
assumed to be constant across the entire range of 
star-forming galaxies. 

The "local" set of constraints is based on in-situ 
measurements of CRs at Earth. These are: 

10 While some of the light absorbed by dust is emitted as far- 
infrared (40 - 120 /im), some is also emitted in near or mid infrared. 
For simplicity, we assume that TIR light is directly proportional 
to the FIR emission . For starburst galaxies, Ltir ~ 1-75-Lfir 
llCalzetti et aI.I2000T l , which we apply to every galaxy for simplicity. 
This corre c tion h as been applied to the reported Lfir/ -kra.rl io in 
lYun et all i 200111 to get our quoted observed -Ltir /^radio • I Belli 
1)2003 ) reports a similar Ltir/^fir ~ 2 for galaxies on the FRC. 

11 Our version of this equation divides our calculated Ltir by 
-^tir/-^fir. ~ 1-75 (Calzetti et al. 2000) to get to the true, ob- 
served FIR emission. See Helou et alj I) 19851 ) for the usual defini- 
tion of q. 



(a) For each electron energy, we calculate the ra- 
tio of the CR positron number density to 
the total number density of positrons and 
electrons at GeV energies. Below GeV en- 
ergies, solar modulation of CRs can af- 
fect the observed CR spectrum. Above 
~ 1 GeV, the observed positron flux 
exceeds the pr edicted flux even in de- 
t ailed models fseelMoskalenko fc Strong 119981 : 
iBeattv et all |2004| lAdriani et al.l 1200a Tut 
see lDelahave et al.ll2009f) . The observed value 
of e + /(e + + e~) i s 0.1 at GeV energ ies (e.g., 
lSchlickersen[2002l: lAdriani et al.ll2009l) . 

We require as a local constraint that 0.05 < 



+ 



e+/(e 
0.0025 g cm' 



-) < 0.2 at 1 



GeV when E ff = 



(b) We also compute the ratio of the proton num- 
ber flux and electron number flux at GeV en- 
ergies. The ratio is observed to be p/e « 
100 at Earth at energies of a few GeV (e.g . , 
IGinzburg fc Ptuskin1ll97^ IS^njTckeiserll2002h . 
This value is also inferred from SN rem- 
nants, which are be lieved to accelerate CRs 
(|Warren et al.H2005T) . 

We require as a local constraint that 50 < 
p/e < 200 at 10 GeV when £ 5 = 
0.0025 g cm- 2 . 

As an alternative way to find the CR proton nor- 
malization, we considered a separate "integrated" 
constraint for the entire Milky Way galaxy using 
an average E g inferred from the Galactic scale ra- 
dius and star formation rate. Our purpose was to 
assess the possibility that the Earth is not in a rep- 
resentative location of the Galaxy; for example, it 
sits in the Local Bubble. 



(a) We calculate the gamma-ray luminosity of the 
Galaxy from tt° decay. We approximate the 
Milky Way as a uniform disk with R — 4 kpc, 
and a surface density of S s = 0.01 g cm -2 de- 
rived from the Schmidt law ( Section 12.21) and 
the Milky Way lum inosity, « 2 x 10 10 L o 
(|Freudenreichl Il998t similar results are ob- 
tained by using th e starlight radiatio n field in 
Strong et al.ll2000L or the S N rate in iFerrierel 
2001h . IStrong et all (|20 0) calculate the to- 



tal 7r gamma-ray luminosity to be L- KQ 
2 x 10 39 ergs s _1 . 

We require as the integrated constraint that 
1 x 10 39 ergs s" 1 < L^ < 4 x 10 39 ergs s _1 
when Eg = 0.01 g cm" 2 . 



Additional checks: As an added check, we have the 
observed CR spectrum at Earth. At hig h energies (7 
1), the observed CR electrons have (e.g.. lLongairlll994l ) 



dh _ cN e (E) 
dE 



47T 



= 0.07 



GeV J 



-3.3 



cm" 2 s _1 sr" 1 GeV" 1 



(15) 
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and the observed CR protons have fe.g.. iMoril Il997t 
iMenn et al]l200(it IAMS Collaboration eTatl 1200^) 



dl p _ cN p (E) 



dE 



47T 



-2.7 



cm 2 s 1 sr 1 GeV 



(16) 

The predicted CR spectrum does not determine whether 
a model was considered formally "successful" , but it was 
used to select among the adequate models for the best 
standard set of parameters. 

Although we did not use it directly, we also calculate 
the spectral slopes at 1.4 GHz from the radio synchrotron 
spectrum as a sanity check. These include the instanta- 
neous spectral slope ot\ 4 as well as the spectral slopes 
to 4.8 GHz (af;!) and 8.4 GHz (af;|). Unless otherwise 
stated, a refers to af;f, the spectral slope from 1.4 GHz 
to 4.8 GHz. Typical values of a are 0.7 — 0.8. As a con- 
straint, a can be very sensitive to minor details in the 
model; we note that a difference of 0.2 in a results in only 
a 60% difference in the specific flux after one decade in 
frequency, and we are mainly concerned with fa ctor of 2 
accuracy in our models. iLisenfeld fc Volkl (|2000D have ar- 
gued that a is decreased by ~ 0.1 by SN remnants within 
galaxies, so our value of a is uncertain at that level. We 
also do not include free-free emission, which can flatten 
the spectral slope, especially in low surface density galax- 
ies. In ULIRGs like Arp 22 0, the observed a is typically 
~ 0.5 (| Clemens et al.l [20081) . but radio emission in these 
galaxies may suffer free- free absorption which flattens the 
spectrum; the un absorbed synchrotr on spectrum a may 
be as high as 0.7 ([Condon et al.| [l991). To some extent, a 
small to moderate difference in a from its observed value 
can be adjusted by altering p, since decreasing p by 0.1 
generally decreases a by 0.05, and p often is not well 
constrained in the considered range 2 — 2.6. Given these 
uncertainties, caveats, and sensitivities in a, and given 
the vast range of galaxies and starbursts we are consid- 
ering, and the simplified parameterizations we are using, 
we do not impose any direct constraint on a. Of course, 
models of individual galaxies should and do account for 
a when they model the radio emission. 

Throughout this work, we assume that the local val- 
ues of the proton normalization and propagation - in 
particular, 6, rj, and / - are the same for both normal 
galaxies and starbursts. We use this assumption for sim- 
plicity, and to keep the number of free parameters rea- 
sonable. In practice, the CR acceleration efficiency and 
the proton-to-electron ratio may change somewhat from 
normal galaxies and starbursts, but we do not consider 
small variations necessary for a basic understanding of 
the FRC. More detailed models of individual systems 
can and do take these changes into account, and we re- 
fer readers to these models if they wish to understand 
starburst galaxies in detail. It is also conceivable that / 
changes dramatically from normal galaxies to starbursts. 
Again, we do not consider this possibility in this paper, 
although we will explore the consequences of very low / 
applying to only starbursts in a future paper. 

3. REVIEW OF PHYSICAL EFFECTS OF PARAMETERS 

To search for models that satisfy the observational con- 
straints listed in Section [2. 3[ our grid of models spanned 
values of a (eq. [T3J), / (eq. [7]), £ and S (eqs. [3] and 2]), 
and p (Section I2.1[) . For a listing of these parameters of 



the model, see Table [TJ As background for interpreting 
our results in SJU we briefly review the effects of these 
quantities on observables. 

3.1. Injection Parameters: 77, 5, and p 

The parameter £ is the normalization of the injected 
primary CR electron spectrum, and with 6 = r;/£, the 
injected CR proton spectrum normalization (eqs. [3] & 
[4]). Changes in £ do not affect the shape of any of the 
equilibrium CR spectra. For fixed 5, larger £ linearly 
increases the CR luminosity and energy density within 
galaxies, and thus — for fixed galaxy parameters — the 
luminosity of CRs in all wavebands, including the radio, 
neutrino, and gamma-ray luminosities. 

An increase in 5 at fixed £ raises the number of secon- 
daries from protons, the ratios e + /(e + + e~) and e~ c /e~, 
and the luminosity from pion decay, L T . 

Of note is the ratio of injected protons to elec- 
trons at relativistic energies. Suppose the electrons 
are injected with a spectrum Q e {E) = C e E~ p and 
the protons are injected with a spectrum Q p (E) = 
C p E~ p . Then, given our normalization condi- 
tions e C R, = X^ me ° C e KE~ p dE and e C R, P = 

J^J" 1 " CpKE-PdE, where K is the kinetic energy 
(see Section 12.11) , it can be shown that 



p-2 



(17) 



The quantity 6 represents the proton to electron ratio 
at high energies (m p c 2 <§C E <C 7maxW e c 2 ) if there were 
no escape, energy losses, or secondary production. Note 
that it is not generally equal to 5, since S is largely depen- 
dent on the shape of the spectrum at low energies. Our 
injection spectra go as E~ p , where E is the total energy: 
the electron spectra stretch down to m e c 2 while the pro- 
ton spectra only extend down to m p c 2 . For steep spectra 
(p > 2), the low-energy particles receive most of the en- 
ergy, so that electrons with E < m p c 2 act as a hidden 
reservoir of energy. 12 This reservoir is unconstrained be- 
cause the observables we use do not constrain the shape 
of the CR spectra at low energies (see Section |2~3"]) . This 
follows from the fact that the FRC is observed at GHz 
frequencies, implying electron energies of order 100 MeV 
to 10 GeV (eq. \T§§. Thus, the actual quantity we con- 
strain is 8. Note that the relationship between C p /C e 
and 6 would be different for another spectrum, such as 
K -p or 7 -p 

The spectral slope p of the injected CRs in part 
controls the final, propagated spectral slope V = 
dlogN (E) / dlogE . The spectral slope, in turn, deter- 
mines how much the secondary particles are diluted. 
Protons at energy E produce secondary electrons and 
positrons of energy E' < E; a steeper primary spectrum 
increases the number of primaries at these lower energies 
compared to the proton energy E. Therefore, a larger 
p (and thus a bigger V for primary electrons) implies a 
smaller e + /(e + + e~) and e~ c /e~. This dilution implies 

12 Conversely, the proton spectrum extends to a maximum en- 
ergy of 7max"ipC 2 , much greater than the maximum energy of the 
electrons; for shallow spectra (p < 2), the reservoir of energy in 
these high energy protons would lower C p /C e at E < 7max>n e c 2 . 



that even in the limit of full proton calorimetry, primary 
electrons may be more important than secondaries. Sim- 
ilarly, the secondary fraction is not a good measure of 
proton calorimetry in itself. For our standard model, 
though, we find that in proton calorimeters, secondary 
electrons and positrons outnumber the primary electrons 
~4-l. 

3.2. Magnetic Field 

The magnetic field strength affects the CR spectra in 
several ways: 

1. It determines the importance of synchrotron cooling 
relative to other radiative and escape losses. The syn- 
chrotron cooling timescale for CR electrons and positrons 
emitting at frequency j^ghz = v/GHz is 

Wnch « 4.5 X 10 7 .Bro'^GHz 2 3*, (18) 

where B\q = B/lOfiG. For normal galaxies, t S ynch is 
comparable to, but somewhat longer than, the inferred 
diffusive escape timescale for the CR electrons producing 
GHz emission in normal galaxies (eq. [5]) . For the ~ mG 
(and larger) fields thought to exist in the densest star- 
bursts, t S y n ch is shorter than even the advection timescale 
(eq.©. 

2. The relative importance of synchrotron also affects 
the propagated equilibrium spectral slope V of electrons 
and positrons; stronger magnetic fields imply steeper fi- 
nal spectra (see Section l3~Tj) . In the limit that cooling 
dominates escape, and that synchrotron is the main form 
of cooling, the equilibrium spectral slope is V = 1 + p. 

3. The magnetic field strength determines the critical 
synchrotron frequency (vc) for electrons and positrons: 

^c«3.3(^) 2 B 10 GHz. (19) 

At a fixed observed frequency (such as 1.4 GHz), a 
stronger magnetic field implies that we see lower energy 
electrons and positrons. 

4. When synchrotron cooling dominates over other 
cooling and escape losses, a stronger magnetic field low- 
ers the equilibrium energy density of CR electrons and 
positrons, because of increased losses. However, in this 
calorimeter limit, each electron and positron has a higher 
luminosity. Therefore, L ra dio approaches a maximum set 
by £, and is not affected by further increases in the mag- 
netic field strength. This effect is the essence of the orig- 
inal calorimeter theory. 

All else being equal in our models of non-calorimetric 
galaxies, larger magnetic fields imply that a larger share 
of injected CR electron power is lost to synchrotron, be- 
cause of the faster synchrotron cooling time. In cases 
when synchrotron does not already dominate, increasing 
the magnetic field strength thus increases Lradio ■ 

Note that the magnetic field in our models is normal- 
ized to the local Solar Circle gas surface density (eq. lT2"|) . 
so that changing a has no effect on local Milky Way con- 
straints discussed in Section [2~3l 

3.3. Effective Density 

The ISM density encountered by CR protons controls 
the production rate of secondary electrons and positrons, 
as well as gamma rays and high-energy neutrinos, from 



inelastic proton-proton collisions. The proton lifetime to 
pion losses is 

^5xl0V(^)", (20) 

from iMannheim k, Schlickeiserl (|1994f) (see also iTorres! 
l2004f ). Higher n e g = f(n) (eq. d) means more secon- 
daries, higher £ ra dio, and higher gamma-ray and neutrino 
luminosities. The secondary electrons and positrons raise 
the ratios e + /(e + + e~) and e~ c /e _ , and they lower the 
equilibrium ratio p/e. Additionally, if the ratio of pri- 
maries to secondaries changes with energy, then the com- 
bined spectral slope V for electrons and positrons can be 
altered, which affects the observed radio spectral slope 
(see Section l3~Tj) . 

The effective ISM density also determines the efficiency 
of bremsstrahlung and ionization losses for CR electrons 
and positrons, with higher densities making these pro- 
cesses more efficient. The bremsstrahlung and ionization 
energy loss timescales are 

tbrcms « 3.7 x 10 7 yr ( , (21) 
\ cm 6 J 

and 

t ion « 2.1 x 10 8 B w ^l 1 yr, (22) 

respectively, where we have again scaled the energy de- 
pendence of iion for CR electrons and positrons emitting 
at GHz frequencies for comparison with t S ynch (eq. 1181) . 
Importantly, energy lost to bremsstrahlung and ioniza- 
tion is not radiated in the radio, so higher / implies lower 
iradio from these processes, all else being equal. 

The energy dependence of these cooling processes also 
flattens the propagated equilibrium electron and positron 
spectra (see the discussion after eq. [T] Section l3Tj) . For 
example, when t sync h = tbrems at some energy and all 
other losses are negligible, then V = p + 1/2 and a — 
p/2 — 1/4. Similarly, when t S ynch = iion and there are no 
other losses, V — p and a — p/2 — 1/2. 

3.4. The Schmidt Law and the Photon Energy Density 

The energy density of photons, and thus the impor- 
tance of IC losses for CR electrons and positrons in star- 
forming galaxies, is set by the slope and normalization of 
the Schmidt law. The IC cooling timescale for CR elec- 
trons and positrons emitting radio synchrotron at fre- 
quency v is 

t IG « 1.8 x 10 8 B\(*v££u£ _ ia yr (23) 

where U p h, -12 = C/ p h/10 -12 ergs cm -3 is the photon 
energy density scaled to that for a typical star-forming 
galax y. For optically-th in galaxies obeying the Schmidt 
law of lKennicutTJ (|l998) and ignoring the CMB, F^/c = 
C^ph,* oc Ssfr <x 4 (eq. [9]). Then, for fixed frequency, 

the IC lifetime therefore scales as t\c °< Sg^ 2-1 ' 4 if the 
magnetic field strength varies as B oc £°. We do not con- 
sider variations on the Schmidt law (e.g.. iBouche et ahl 
2007|), but their effects can be inferred from equation 1231 
if Ssfr has a steeper increase with S g , then tic wm fall 
more rapidly with surface density. In practice, the CMB 
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will make IC losses more efficient in the lowest density 
galaxies, and any FIR opacity (Section 12.2.4ft will make 
them more efficient in high-density starbursts. 

Just as bremsstrahlung and ionization losses can re- 
duce the share of energy left for synchrotron radiation, 
a greater photon energy density and IC power decreases 
^radio- Unlike bremsstrahlung and ionization, though, 
IC losses produce a steep spectrum. In the limit where 
they dominate other losses and escape, V = I + p, as in 
the case of pure synchrotron cooling, and a = p/2. 

4. RESULTS 
4.1. Standard Model 

We adopt p = 2.3, / = 1.5, a = 0.7, 8 = 48 (5 = 5.0), 
and £ = 0.023 as our fiducial model. This model re- 
produces the FRC, as seen in Figure [1] (solid line). In 
this particular model, we require £ = 0.023 to match the 
normalization of the FRC. The ratio of FIR to 1.4 GHz 
luminosities varies by only 1.7 over the entire range of 
E g , and shows no obvious trend. However, the scatter 
appears to be concentrated at the low-E g end of the FRC, 
with iTm/iradio varying by less than 12% in the star- 
bursts in this model. 

The standard model also satisfies both local and inte- 
grated constraints on the proton normalization, as well 
as the observed CR spectrum. Our positron ratio at 1 
GeV, e + /(e + + e~) = 0.10, and proton-to-electron ratio 
at 10 GeV, p/e = 82, are good matches to the observed 
values. The Milky Way 7-ray luminosity in this model, 
2.0 x 10 39 er g s s" 1 , is also a good match to the value of 
IStrong et all (|2000T) . 

The predicted proton spectrum at Earth in this model 
is 91%-118% of its observed value at 1, 10, and 100 GeV, 
implying that p is well matched to the Galactic CR spec- 
trum. Similarly, the CR electron flux at Earth is 122% 
of its observed value at 10 GeV. The least satisfactory 
aspect of this model is the predicted spectral slope a for 
Milky Way-type galaxies (a ps 0.9 — 1.0), which is some- 
what too high. Our results are, however, reasonable for 
starbursts (a ps 0.5 — 0.7; see Figure [2]). 

We emphasize that the parameters of our standard 
model are adequate for all star-forming galaxies on the 
FRC. We discuss the many competing effects that yield 
the FRC in Section [Ol 

4.2. Degeneracy in the Standard Model 

Our local set of constraints ( Section |2. 3ft narrow down 
the allowed parameter space considerably. The models 
that survive have 6 ps 34 — 100 and a = 0.6 — 0.7. To 
get the correct normalization of the FRC, we must set 
0.019 < £ < 0.027 when p = 2.3, so that 0.097 < rj < 
0.22. Flatter injection spectra generally have lower £ 
(down to 0.006 for p = 2.0) and higher r\ (reaching 0.28 
for p = 2.1), while steeper injection spectra generally 
have higher £ (up to 0.18 for p — 2.6) and lower r\ (as 
low as 0.09, which occurs when p = 2.4). However, ot\\ 
is somewhat high when p > 2.2 (0.9 — 1 predicted com- 
pared to 0.7 — 0.9 observed) for normal galaxies, but is 
close to observed values for p < 2.2 (0.8 — 0.9 predicted). 
The spectral slope is sufficiently low (af'f ~ 0.5 — 0.7 
predicted) for starbursts. 

The integrated Milky Way 7-ray luminosity from ir° 
decay provides similar, but somewhat weaker constraints, 



favoring lower 6. At low p = 2.0, models with 1 < 5 < 
50 are selected by the 7-ray luminosity. Higher p models 
continue to work so long as S decreases, because the nor- 
malization depends on the spectrum at very low energies 
as discussed in Section 13.11 We can take the normal- 
ization into account by comparing S (eq. 1171) . and we 
find that the allowed S (10 < S < 100) slowly increases 
with p (from roughly 25 at p — 2.1 to 91 at p = 2.6 
when / = 1.0) and decreases with / (from roughly 60 at 
p = 2.2, / = 1.0 to 50 at p = 2.2, / = 2.0). Even p = 2.6 
models predict an FRC and the correct tt° luminosity of 
the Milky Way; we would need to take into account either 
the observed a in normal galaxies or the local observed 
CR spectral slope to further constrain p in our standard 
model. For example, p = 2.6 works when 6 ss 90 (5 ps 1), 
and £ ps 0.17. As discussed in Section ETT1 a higher p di- 
lutes secondaries and lowers the fraction of electrons that 
are secondaries. Therefore, the secondaries contribute a 
smaller fraction of the radio luminosity and are less likely 
to break the FRC as galaxies become proton calorime- 
ters at high density. The FRC can then tolerate a higher 
secondary production rate (and ultimately higher 6) for 
high p. A higher £ is needed, though, since more of the 
energy goes into unobserved low energy electrons. These 
high p solutions also produce steep synchrotron spectra 
with a ps 1.1 — 1.2 in normal galaxies, and can be con- 
strained by the spectral slopes. 

More broadly, we also consider variations on our usual 
parametrization, such as lifetimes including advection, 
different scale heights, and FIR optical depths. Many 
of these variations are inconsistent with the constraints 
in Section 12.31 However, those that do satisfy the con- 
straints had similar values for £, 5, p, and a as the fidu- 
cial model. We describe these variants in detail in Ap- 
pendix [X] 

4.3. General Features of the Particle Spectra 

We show typical predicted CR spectra in Figure [3] for 
Yj g = 0.001 g cm~ 2 and 10 g cm~ 2 , the lowest and highest 
surface densities we consider. 

In low surface density galaxies, the protons with 7> 1 
have a power-law spectrum with V p about 0.5 greater 
than the injected spectral index p. The increased steep- 
ness co mes from faster diffusive esc ape at higher energies 
(eq. iGinzburg fc Ptuskinl 1 197a) . At lower energies, 
the CR proton spectrum flatten s due to ionization losses, 
which are constant with energy (|Schlickeiserll2002l : iTorresI 
|2004| ). High surface density galaxies have harder proton 
spectra with V p = p at high energies, since pion losses 
overwhelm escape and are roughly energy independent. 

In low surface density galaxies, the primary electrons 
behave similarly to protons. For most low energies, they 
have a power-law spectrum with V e ps 0.5 + p, caused b y 
diffusive escape losses (e.g.. IGinzburg fc Ptuskinl 1 1976| ). 
However, synchrotron and IC losses steepen the spectrum 
at high energies. Bremsstrahlung and ionization flatten 
the spectrum at l ower energies (E < 1 GeV) (compar e 
Figure 12 see also iThompson et all 12005 ICondonlll992[ ). 
In high surface density galaxies, diffusive losses are negli- 
gible compared to synchrotron and I C losses, thus forcing 
Vp pa 1 +p at higher energies (e.g.. IGinzburg fc Ptuskinl 
fl976T) . 

The spectra of secondary electrons and positrons show 
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Fig. 1. — The non-thermal FRC, as reproduced in our standard model (p = 2.3, / = 1.5, a = 0.7, S = 5, £ = 0.023). While low CR 
escape times and low UV optical depth on their own would break the correlation at low surface densities, the two effects cancel each other 
out, creating a largely linear FRC. 



additional features with respect to the primary elec- 
tron spectrum. The secondary (pion-produced) spec- 
trum is flatter than the primary spectrum at low en- 
ergies, because the production cross sections decreases 
near the pion production threshold (see Figure [5J 
IStrong & MoskalenkdH99llTorred l2004). At high ener- 
gies, the pion electrons and positrons are injected with a 
spectrum proportional to the steady-state proton spec- 
trum. This means that in low surface density galax- 
ies, the pion electrons and positrons have a steeper 
spectrum than the primar y electrons at high energies 
(|Ginzburg fc Ptuskinll 19761) , while in high surface density 
galaxies, the secondary and primary spectral slopes are 



the same at high energies (as seen in Figure [5]). Further- 
more, there are always more secondary pion positrons 
than secondary pion electrons (ultimately due to charge 
conservation). Knock-off electrons become increasingly 
important at ver y low energie s and dominate as 7 ap- 
proaches one (see lTorresll200l . 

We also show radiation spectra in Figures |5] and [7] for 
synchrotron, ir° 7-rays, relativistic bremsstrahlung, and 
IC emission (see Section 1531 for the assumptions used to 
estimate the IC emission). The radio, bremsstrahlung, 
and IC emission generally steepen with increasing fre- 
quency (compare with Figure HJ see also iLisenfeld et al.l 
ll996atlThompson et al.ll2006l) . whereas pion 7-rays peak 
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Fig. 2. — Spectral slope as a function of the gas surface 
density. In this plot, 01.4 is the instantaneous spectral slope, 
dlog Fjs/dlog v, at 1.4 GHz. Elsewhere in the paper, a is the 
observable ct 4 -|. The parameters have their fiducial values. For 
our standard p = 2.3, strong cooling by synchrotron and IC alone 
would imply that a = 1.15. Instead the spectral indices are signif- 
icantly flatter, especially in the strong cooling calorimeter limit at 
high S 9 , as a result of ionization and bremsstrahlung losses. 

at a few hundred MeV. Although we do not calculate it 
here, the overall high-energy neutrino emi ssion is compa - 
rable to the 7-ray emission from 7r° decay (|Steckerlll979t ). 

5. DISCUSSION 
5.1. Is Calorimetry Correct? 

We show the effects of forcing electron and UV 
calorimetry to hold in Figure [T] for our standard model 
(cf. Section I4.1j) . It is clear that most of the energy 
in 1.4 GHz electrons is lost radiatively in galaxies with 
E 9 > 0.01 g cm~ 2 : calorimetry holds in high- but not 
lo w-density galaxies (this behavior was first described 
in IChi fe; W olfcnd alel [l990l and was also predicted by 
iLisenfeld et al.lll996a|) . At lower surface densities, elec- 
tron calorimetry begins to fail (decreasing the radio lu- 
minosity), but the effect of this on the FRC is largely 
mitigated by the decreasing optical thickness to UV pho- 
tons (decreasing the FIR lumi nosity) . T his conspiracy 
saves the FRC, as discussed by Bell (2003), but only ap- 
plies in our standard model over one decade in E ff . In 
our standard model, electron escape at low E g eventually 
becomes the stronger effect, so that low-density galaxies 
would be radio dim with respect to the FRC. A high value 
of q is in fact ob served as a nonl inearity in the FRC at 
low l uminosities (|Yun et al.ll2001l though [Beswick et al.l 
120081 find the opposite). Unfortunately, studies of the 
low luminosity FRC are com plicated by the pre sence of 
thermal radio emission (e.g.. (Hughes et al.ll2006l for the 
Large Magellanic Cloud), which also correlates with FIR 
light and overwhelms the nonthermal synchrotron emis- 
sion considered here in the lowest density galaxies. 

While the standard models predict that electron 
calorimetry holds in the inner Milky Way, there are phys- 
ically motivated variants (Appendix which predict 
that electron calorimetry fails for normal galaxies and 



the weakest starbursts. In particular, the "strong wind" 
variants predict a non-calorim etric inner Mil k y Way , be- 
cause of the wind inferred by lEverett et al.l ( 20081 ) (see 
Appendix IA.2[) . 

Although the transition to calorimetry is model de- 
pendent, it seems unavoidable that extreme starbursts 
like Arp 220 are electron calorimeters. We can derive 
the speed v esc at which CRs would have to stream out 
of galaxies for electron calorimetry to fail, according to 
the cooling rates in our standard model. We can also 
compare these numbers to standard CR confinement the- 
ory, where CRs are limited to propagate at the Alfven 
speed (va = B/yJAirp) by a streaming instability in the 
ionized ISM (jKulsrud fc Pearcd [19691 ) . We can also in- 
vert the problem and determine the magnetic field with 
a high enough Alfven speed 13 for CRs to stream out 
of the galaxy in one cooling time (B 2 SC = 47T/w 2 sc ), 
as well the diffusion constant (D csc = h 2 /t coo \) needed 
to diffuse out of the system in one cooling time. The 
environmental conditions needed to allow CRs to es- 
cape before cooling significantly are reasonable for weak 
S g < 0.1 g cm~ 2 starbursts. We find that v csc = 
620 km s _1 , and winds of several hundred kilometers per 
second are in fact observed in starbursts. Similarly, we 
calculate D csc = 1.9 x 10 28 cm 2 s _1 , and diffusion con- 
stants of ord er 10 28 cm 2 s -1 are inferred for starburst 
galaxies (e.g.. iDahlem et aT1ll995| ). However, if the CRs 
stream through mean density ISM, then B csc ss 3 mG 
is higher than the equipartition magnetic field strength 

8tt 2 GE 2 « 300 fiG for E g 



B, 



0.1 g cm , 

so that CR escape would have to be super- Alfvenic. 
For higher E 9 , though, escape would require extreme 
wind speeds (8000 km s _1 when E g = 



0.4c when E„ = 10 g cm -2 

-2 



1 g cm 2 and 



extremely 

2 



for E s = 1 g cm 



120, 000 km s" - „„. „ 
high diffusion rates (2.5 x 10 29 g cm 

and 3.6 x 10 30 g cm -2 for E g = 10 g cm -2 ), or extremely 
strong magnetic fields (B csc = 0.11 G in the E g = 
1 g cm -2 case and B csc = 5 G in the S g = 10 g cm~ 2 
case) that are unreasonable. We therefore conclude that 
electron calorimetry must hold in dense starbursts. 

We can similarly ask whether galaxies are proton 
calorimeters. The low pion luminosity of the Galaxy and 
the secondary positron fraction at Earth imply that nor- 
mal galaxies like the Milky Way are not proton calorime- 
ters. We have estimated the proton calorimetry frac- 
tion F ca x in our models by adding the emissivity in pion 
products 14 to the emissivity in CR protons with energy 
greater than 1.22 GeV, the pion production threshold 
energy. However, when we do this we find that even ex- 
plictly proton calorimctric models with no diffusive or ad- 
vective escape have F ca ,\ sa 0.5. This appears to be caused 
by an inconsis tency between the pionic l ifetim e we use 
(cqn |2T)|) from Ma nnheim fc Schlickeiserl (|1994[ ) and the 
GALPROP cross sections: if we add up all of the energy 
in all of the pionic products of a CR proton of energy ~ 
GeV, the effective energy loss rate is several times smaller 

13 This estimate assumes that CRs are streaming through mate- 
rial with the mean ISM density. In a lower density phase, va wm 
be larger and B csc will be smaller. 

14 Since we do not calculate the neutrino spectrum, we simply 



assume that Q v 
high energies. 



Q 7 , which is a reasonable approximation at 
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Fig. 3. — The predicted kinetic energy spectra of cosmic rays in a low density galaxy (left, E 3 = 0.001 g cm -2 ) and a high density 
starburst (right, S s = 10 g cm -2 ) for our standard model (p = 2.3, / = 1.5, a = 0.7, 5 = 48, § = 0.023). We mark the kinetic energy 
where electrons and positrons emit synchrotron radiation at v « 1.4 GHz. The cutoff in the lepton spectra at ~ 500 GeV is caused by our 
use of 7 m ax = 10 6 ; see Section 14.31 for further discussion. 




Fig. 4. — Left: cooling times (for electrons and positrons) when vq = 1.4 GHz. Right: emissivity (energy lost per volume per unit time) 
for each process from protons, electrons, and positrons, integrated over energy. Pion losses include all of the energy going into secondary 
production as well as 7-rays and neutrinos. 
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Fig. 5. — Left: the abundance of positrons as a function of energy. Right: the abundance of secondary electrons compared to all electrons, 
as a function of energy. The saturation of the positron and secondary elect ron ratios at high S a , and the lack of energy dependence for 
E > 1 GeV, is a sign of proton calorimetry (see Section 14.31 and Section 15.11 1. We mark with filled symbols the energies where the critical 
synchrotron frequency vc is 1.4 GHz. 
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Fig. 7. — The spectra of bremsstrahlung 7-rays (left) and estimated IC (right) predicted by our standard model (p = 2.3, / = 1.5, 
= 0.7, <5 = 48, £ = 0.023). We discuss our assumptions for the IC emission in Section 15.51 

than implied by IMannheim fc Schlickeiserl (119 94) 15 . Wc 
also note that the IMannheim fc Schlickeiserj (1994) pio- 
nic lifetime is twice as shor t at ~ GeV energies as the pio- 
nic lifetime in Schlickciscr (2002). Finally, this approach 
ignores ionization losses, which do not create secondaries 
but will prevent lower energy protons from escaping. 

To account for this discrepancy in the energetics, we 
normalize our estimate of F ca \ so that an explicitly pro- 
ton calorimetric model of the same CR injection rate, £ s , 
p, and / has F ca \ = 1. We then see in Figure[5]that dense 
starbursts with E 9 m 10 g cm" 2 all are proton calorime- 
ters with F ca \ R3 1 for several variants (Appendix |^J) . 
As with electron calorimetry, proton calorimetry some- 
times breaks down for the E s = 0.1 g cm" 2 weak star- 
bursts, because the time to cross the 100 pc starburst 
scale height is short. However, when E 9 > 1 g cm" 2 , 
proton calorimetry holds in our models; a model with 
winds and strong diffusive losses has proton calorimetry 
breaking down at S g = 1 g cm" 2 (F cai \ — 0.45). 

As with the electrons, we can derive the speed that 
the CRs would need to stream out of a starburst for 
proton calorimetry to fail, which is v csc — h/t n — 
1900 km s (£ 9 /g cm -2 )/. While v csc is easily at- 
tained by winds in starbursts with E s = 0.1 g cm" 2 , 
only the fastest winds are capable of breaking pro- 

1 g cm -2 starbursts. Dif- 
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Fig. 8. — The estimated proton calorimetry fraction F ca i in our 
models for several variants. The values are normalized so that 
models with no escape have -F ca i = 1. Normal galaxies are not 
proton calorimeters, while ULIRGs with S 9 « 10 g cm -2 are in all 
of our variants (note the convergence of all models to ~ 1 at high 
S 9 ). In all variants, proton calorimetry holds for starbursts with 



2 . Variants shown are our standard model (solid 
1.5, 6 = 48, a = 0.7); B oc p a with 300 km s" 1 



_ 1 g cm 
p = 2.3, / = 

wind in starbursts (dash; p = 2.2, / = 1.0, 8 
strong winds of 175 km s" 1 in S 9 = 0.01 g cm 
for starbursts (long dash dot; p = 2.2, / = 2.0, S = 45, a = 0.5) 
constant D z , B oc p a , and winds of 300 km s _1 in starbursts (long 
dash; p = 2.2, / = 1.5, 8 = 34, a = 0.6); and fast diffusive escape 



= 67, a = 0.5) 
~ 2 , 600 km s _] 



with B oc 
a = 0.6). 



S? and no winds (dotted; p = 2.2, / = 2.0, 5 



45, 



ton calorimetry in E s 
fusive escape limited to the mean Alfven speed of the 
starburst would require strong magnetic fields (i? C sc = 
27 mG(/i/100pc)- 1 / 2 (E 9 /g cm" 2 ) 3 / 2 /) with energy den- 
sities greater than the midplane gas pressure in star- 
bursts to break proton calorimetry. We therefore con- 
clude that proton calorimetry is difficult to avoid in / w 1 
starbursts with E s > 1 g cm" 2 . 

15 As far as we are aware, this discrepancy has not been discussed 
in the literature. 
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Fig. 9. — The high-E s conspiracy in our standard model (p = 2.3, 
/ = 1.5, a = 0.7, 5 = 48, f = 0.023). The simple calorimeter model 
has perfect UV calorimetry and electron calorimetry, with only 
synchrotron cooling and no secondaries. Non-synchrotron cooling 
and secondaries alone each create a broken FIR-radio correlation, 
but conspire to make it linear at high density. 

5.2. What Causes the FIR-Radio Correlation? 
5.2.1. Calorimetry and the vc Effect 

Calorimetry provides a simple way to explain the FRC. 
We find that both electron and UV calorimetry hold 
for starbursts, and possibly the inner regions of normal 
galaxies, depending on the variant on our underlying 
model (Appendix IA"|) . Calorimetry therefore serves as 
the foundation of our explanation for the FRC. Other ef- 
fects alter the radio luminosity, both at low density and 
high density, but by a factor of ~ 10, compared to the 
dynamic range of 10 4 in E ff . At the order-of- magnitude 
level, calorimetry can be said to cause the FRC, and 
other effects are relatively moderate corrections. 

However, in more detail, we find that iTm/^radio is 
not in fact flat even in the simple calorimeter model, 
with no escape, non-synchrotron cooling, or secondaries 
(the light dotted line in Figure H]). Instead, iTm/^radio 
decreases by a factor of 2.6 as E s increases, because 1.4 
GHz observations probe lower CR electron energies as 
the magnetic field strength increases. We call this de- 
crease in ixm/iradio with E s the "vc effect". In gen- 
eral, the effect becomes more significant as p increases 
past 2.0, because the electron spectrum becomes steeper. 
It can be shown that in this simplest calorimeter limit, 

iTIR.AradioCX.BP/ 2 - 1 . 

5.2.2. High-Tig Conspiracy 

The radio luminosity in high-density galaxies is altered 
from the calorimetric luminosity mainly by two mecha- 
nisms, non-synchrotron cooling and the appearance of 
secondary electrons and positrons. We illustrate these 
effects in Figure^ 

In normal galaxies, synchrotron cooling dominates the 
energy losses, though bremsstrahlung and IC off the 
CMB can be competitive within a factor of a few or 
less. However, in starbursts, energy loss is mainly by 



bremsstrahlung and ionization. This decreases the pro- 
portion of energy lost that goes into radio. The energy 
diverted to bremsstrahlung and ionization therefore in- 
creases iTm/^radio by a factor of up to ~ 20 in star- 
bursts compared to normal galaxies (compare the dotted 
and long-dashed lines in Figure |9]) . 

Secondary electrons and positrons themselves radi- 
ate in the radio. In the starbursts, which are proton 
calorimeters, there are several times more secondaries 
than primary electrons, while in normal galaxies, the sec- 
ondary contribution is small. Secondaries increase the 
radio emission by a factor of ~ 4 in starbursts compared 
to normal galaxies (compare the dotted and short-dashed 
lines in Figure [5]). 

These effects each on their own alter the calorimetric 
radio luminosity by up to an order of magnitude. Since 
both are density dependent, they both become impor- 
tant in starbursts. However, combined with the vc ef- 
fect (Section r5.2.1|) in the simple calorimeter model, they 
largely cancel each other out to maintain a linear FRC. 
The exact magnitudes of these effects are model depen- 
dent, but they are always important and the direction 
each works in is the same in every case. It is possi- 
ble that relaxing the assumptions of our approach, such 
as including time dependence or spatial variation, could 
avoid the severe non-synchrotron losses and secondary 
electrons and positrons giving rise to this particular high- 
ly conspiracy. However, any new effects would have to 
be tuned to avoid the processes we already include while 
still reproducing the FRC, trading one conspiracy for an- 
other. 

There are two other effects that appear in our vari- 
ants (Appendix |A"| . but not our standard model, which 
can change the FRC. First, if the magnetic field is as- 
sumed to depend on density instead of surface density 
(Section lA.lj) . the magnetic fields will be much stronger 
in the starbursts for the B <x p a case, since the starbursts 
are more compact. This will make synchrotron cooling 
dominant again, upsetting the high-E g conspiracy. This 
effect can be compensated by winds and a weak magnetic 
field dependence on p (low a). Second, if the FIR opti- 
cal depth is significant (Section I A. 4[) . the photon energy 
density inside the galaxy is greater by a factor of ~ tfir 
than inferred from the photon flux alone. While typical 
FIR opacities are small, the optical depth is appreciable 
in dense starbursts (1 cm 2 g _1 < kfir ^5 10 cm 2 g _1 ). 
This increases IC losses dramatically at the high densi- 
ties, decreasing the radio luminosity. Models with large 
FIR optical depths have trouble reproducing the FRC 
(Section rO|) . 

5.2.3. Low-'Eg Conspiracy 

The radio luminosity in low density normal galaxies 
is modified by a different pair of opposing mechanisms, 
the failure of electron calorimetry and the failure of UV 
calorimetry. This conspiracy is illustrated in Figure [TJ for 
our standard model. 

Normal galaxies are not generally electron calorimeters 
- both diffusive and advective escape can operate faster 
than cooling. In weak starbursts (£ s < 0.1 g cm -2 ), 
escape can be competitive with cooling processes, but 
not in stronger starbursts. Escape therefore decreases 
the radio emission in normal galaxies compared to the 
calorimetric expectation. 
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However, normal galaxies are generally not UV 
calorimeters either; a substantial fraction of the UV light 
emitted by star formation can e scape withou t bein g re- 
processed into FIR light ("e.g.. IXu fc Buatl Il995t iBeli 



20031: iBuat et aLll2005t Martin et alJ l2005: Po pescu et al 
2005T ) . Therefore, normal galaxies also have a lower FIR 
luminosity compared to the calorimetric expectation. 

As can be seen in Figure [T] each of these effects alters 
the FRC by a factor of — 4 for S g = 0.001 g cm" 2 in our 
standard model. Since they work in opposite directions, 
the resulting iTm/iradio nonetheless remai ns the sam e 
as the calorimetric prediction (as suggested by Bell 2003). 

5.2.4. The Intermediate Case 

The boundary between these two conspiracies occurs 
when E ff = 0.1 g cm -2 . In some variants (Appendix IA"|) . 
factors from both surface density regimes must be tuned 
to maintain the FRC at this surface density: escape time, 
secondaries, non-synchrotron cooling, and magnetic field 
strength all have an effect on the radio luminosity. How- 
ever, these starbursts are unavoidably opaque to UV 
light, so the full low-£ s conspiracy cannot work for these 
galaxies. This becomes a problem when CR escape is 
quic k, su ch as when strong winds are present (see Sec- 
tion causing these galaxies to be radio-dim. Since 
the conspiracies begin to break down for the weakest star- 
bursts, the transition from normal galaxies to starbursts 
may prove important in testing models of the FRC. 

5.2.5. Summary 

The many factors described above conspire to produce 
the FRC, both in low-density non-calorimetric galaxies 
and high-density calorimetric starbursts. The traditional 
distinction between calorimeter and conspiracy explana- 
tions of the FRC is not clear cut in our models. We find 
that the FRC requires both calorimetry and conspiracy. 

5.3. The FIR-Radio Correlation at Other Frequencies 

We have mainly considered the well-studied FRC at 
1.4 GHz. However, the F RC is als o known to exist at 
150 MHz (ICox et al.lll988h . 4.8 GHz (Ide Jong et al lfl985t 
iWunderlich et all 11987ft . and 10.55 GHz (iNiklasI 119971) . 
The correlation holds for both normal galaxies and star- 
bursts at these frequencies, and remains even after ther- 
mal radio emission is subtracted. We show the predicted 
ratios of FIR to synchrotron radio fluxes in Figure [TU1 

Our standard model predicts increased nonlinearity at 
other frequencies for a set of galaxies that span from nor- 
mal galaxies to starbursts. While iTm/^radio varies by 
only 1.7 at 1.4 GHz over the full range in E 9 , it varies 
by 2.3 at 500 MHz, and a factor of 5.4 at 100 MHz. At 
higher frequencies, the situation is similar, though the 
linear FRC is somewhat better preserved: £TiR/£radk> 
varies by a factor of 2.3 at 4.8 GHz, 2.8 at 8.4 GHz, 
and 3.6 at 22.5 GHz. As can be seen in Figure [TUl at 
low frequencies the FRC is predicted to tilt to the FIR 
with increasing S g , while at high frequencies the cor- 
relation is predicted to tilt to the radio in starbursts. 
Our fiducial model with winds and B cx p a (Section fA.2|) 
predicts a similar increase in scatter at other frequen- 
cies (£ T iR/£ r adio varies by 2.1 at 500 MHz; L Tm / L radio 
varies by 2.4 at 4.8 GHz). 

Our models also predict that the normalization of the 
FRC should change with the observed frequency. In 
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Fig. 10. — The FIR-radio correlation at other frequencies. 
Galaxies toward the top of the plot are radio-dim for their FIR 
emission, while galaxies toward the bottom of the plot are radio- 
bright. The shown frequencies are 100 MHz (light gray, solid), 
500 MHz (light gray, long-dashed), 1.4 GHz (medium gray, long 
dashed/short dashed), 4.8 GHz (black, dash dot), 8.4 GHz (black, 
short dash), and 22.5 GHz (black, dotted). All lines are for our 
standard model (p = 2.3, / = 1.5, a = 0.7, 5 = 48, £ = 0.023). 

general, £Tm/£radio decreases with increasing frequency. 
This effect is stronger for the starburst galaxies, where 
the nonlinearities in the predicted FRCs appear. The 
radio-brightness at high frequencies is a direct conse- 
quence of the strong bremsstrahlung and ionization cool- 
ing in our models: synchrotron losses are more efficient 
relative to bremsstrahlung and ionization at higher ener- 
gies, so that more energy goes into radio emission. Only 
when a„ reaches 1 does the radio emission begin to de- 
crease with frequency. 16 

Direct comparison between our models and observa- 
tions can be difficult, because the FRC is usually con- 
sidered in terms of luminosity rather than S g and be- 
cause the FRC is often fit as a nonlinear function. We 
can nonetheless make some qualitative comparisons be- 
tween observations and our models. The observed 151 
MHz correlation appears to be nonlinear, with lumi- 
nous galaxies being bright er in the radio t han would be 
predicted from the FIR (|Cox et al.l 119881 find £fir oc 
L°- 87±0 04 ). Our models predict the opposite effect if 
Ltir increases mono tonically w i th £„ , with iTm/^radio 
increasing with E g . iFitt et al.l (|1988) attribute the ob- 
served non-linearity at these frequencies to the FIR emis- 
sion of old stars, and infer a linear FRC when they re- 
move this effect. It is also worth noting that Arp 220 is 
radio dim at 15 1 MHz, though this may be due to free- 
free a bsorption (jSopp fc Alexander! 119891 : ICondon et al.l 
119911 ). At 4.8 GHz, the FRC is known t o be tight (0.2 
dex dispersion) and approximately linear (jde Jong et al.l 

16 This is also a generic prediction if there are loss processes 
that dominate synchrotron at low energies. For example, galax- 
ies are radio dim at low frequencies if they have strong diffusive 
losses (tdiff oc i? -1 / 2 ) or winds (t w i n d constant with E). A large 
^Tirt/^radio also arises if there is radio absorption at low frequen- 
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Fig. 11. — The instantaneous spectral slope of the (nonthermal) 
synchrotron emission as a function of frequency. In this plot, ct„ is 
the instantaneous spectral slope, dlog F v /dlog v. Elsewhere in the 
paper, a is the observable af'f. 

fl985t iWunderlich et al.lfl98l . though our models pre- 
dict that starbursts should be radio bright compared 
to their FIR fluxes at these frequencies. At 10.5 GHz, 
most of the radio emissio n is thermal and not from 
synchrotron. iNiklasI (|1997l) estimates the contribution 
from synchrotron alone and finds a nonlinear dependence 
L v cx Lpj 2 ^ 09 , so that the FRC tilts towards stronger 
radio emission at higher luminosities. Assuming that 
Ltir increases with E s , we find a qualitat i vely si milar 
behavior in our models. However, INiklasI (]l997|) also 
finds a non-linear FRC at 1.4 GHz, with no dependence 
on frequency for the slope of the FRC, in contrast to 
lYun et al. (2001) who find a linear correlation (except at 
low luminosities) but only consider the FRC at 1.4 GHz. 

At least two effects we do not include would complicate 
our predictions. At low frequencies, free-free absorption 
may significantly lowe r the radio flu x beyo nd what we 
predict in starbursts. ICondon et al.l (|1991|) argue that 
free-free absorption is important even at GHz frequen- 
cies in starbursts like Arp 220, and it becomes more 
effective at low frequency. This effect would make the 
low frequency FRC even more nonlinear than we predict. 
Thermal emi ssion becomes significant at high frequencies 
(|Niklaslll997l estimates that ~ 30% of the radio emission 
is thermal at 10.5 GHz). While the thermal contribution 
can be estimated and subtracted off, at very high fre- 
quencies it may so overwhelm the synchrotron radiation 
that studying the correlation between FIR and nonther- 
mal radio becomes impossible. 17 

5.4. The Spectral Slope a 

In the Milky Way, the observed spectral slope a in- 
creases (the spectrum steepens) with frequency. At low 
frequencies (< 100 MHz), the spectral slope is only 

17 A linear correlation between thermal r adio emission and ra- 
dio e mission is predicted and observed (e.g.. ICon don 1992; INiklasI 
1997), though it provides no information on the cosmic rays or 
magnetic fields in a galaxy and is beyond the scope of this paper. 



w 0.4-0.5 fe.g.. lAndrewlll966tlRogers fc Bowmanll2008ft 
but a reaches « 0.75—0.8 a t GHz frequencies and reaches 
ps 0.8-0.9 at several GHz (|Websterill97llPlatania et al.l 
[19981 H00l before free-f ree emission flattens the spec- 
trum (e.g., iKogut et al.ll2009T ). though there are vari- 
atio ns with direction an d Galactic latitude (for exam- 
ple, EichXRiicE EHI). In fact, our models predict 
a steepening with frequency, though in our standard 
model a. v is higher than observed at all frequencies for 
Eg = 0.0025 g cm" 2 : it is 0.79 at 100 MHz, 0.91 at 1 
GHz, and 1.00 at 10 GHz (Figure HU see also Figures [6] 
and [2]). We note that a can be decreased by adjusting 
p; a value of p = 2.1 can decrease a by 0.1. We also 
note that we used an escape time that increased as en- 
ergy decreased; if the escape ti me is constant or even 
decreasing at low energ ies (e.g.. lEngehnann et al.lll990t 
I Webber fc Higbie][200l . then our low frequency a will 
also decrease, and be more in line with observations. The 
predicted a v for E ff = 0.01 g cm -2 are somewhat better 
at low frequencies: 0.73 at 100 MHz, 0.93 at 1 GHz, and 
1.07 at 10 GHz. Models with p w 2.0 - 2.1 do a better 
job of matching the observed spectral slopes of the Milky 
Way. 

As can be seen by the cooling and escape times in 
Figure [H our standard model implies that escape, syn- 
chrotron, and bremsstrahlung all can shape the spectrum 
in normal galaxies. Escape dominates at low surface den- 
sities, while all three are comparable for the inner regions 
of galaxies (E g 0.01 g cm -2 ). Our greatest problem 
with a for normal galaxies is that it is predicted to in- 
crease slightly with E ff (Figure^. This would imply that 
the inner regions of spirals would have steeper spectra 
than the outer r egions, when in fact the opposite effect 
is observed (e.g.. iMurgia et ai1l2005[ ). 

The reason for the steepening is that escape becomes 
less effective as the galaxies become denser, so that cool- 
ing prevails. In normal galaxies, synchrotron dominates 
bremsstrahlung by a factor of a few (and ionization by an 
order of magnitude) , and the ratio of the synchrotron to 
bremsstrahlung cooling times is only weakly dependent 
on density (compare the short-dashed synchrotron and 
the long-dashed bremsstrahlung lines in Figure 0]). For 
a constant scale height so that B cx (n) a , we have from 
equations [TH] and [3T] that £ S ynchAbrems oc (n)B~ 3 / 2 cx 
(n)i~ 3 °/ 2 , whjch j s essentially constant for a — 0.7 and 
slowly changing for a = 0.5. In contrast, from equa- 
tions [5j HH and [19j we find that at fixed frequency 
t sy nch/*diff oc 5~ 7 / 4 cx (n)~ 7a / 4 , roughly inversely pro- 
portional to (n). This implies that synchrotron losses 
become much more effective than escape as density in- 
creases, but the bremsstrahlung losses remain a factor of 
a few less important than synchrotron losses. Therefore, 
as normal galaxies become calorimetric, their radio spec- 
tra will become steep in our model, since bremsstrahlung 
is only important enough to flatten the spectrum from 
its pure synchrotron-cooled limit of a v = p/2 w 1.1 to 
a v « 0.9 - 0.95. 

This problem remains for all of the variants (Ap- 
pendix |Aj that satisfy local or integrate d co nstraints, 
except in the strong wind variant (Section IA.2[) in which 
we include adve ctive escape that wo uld result from the 
wind inferred by lEverett et al. (2 008), and the fast diffu- 
sive escape variant (Section IA.6[) . If normal galaxies all 
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Fig. 12. — The predicted synchrotron radio spectra of the star- 
bursts in M82, NGC 253, and Arp 220, compared with Allen Tele- 
scope Array observations from Williams & Bower ( 2010). Our fidu- 
cial model is in black (solid - using Schmidt law and estimated 
starburst volume; dashed - scaled to 1.4 GHz flux) and our fidu- 
cial model with winds and B oc p a is in gray (solid - using Schmidt 
law and estimated starburst volume; dotted - scaled to 1.4 GHz 
flux). We do not include free-free absorption or thermal emission. 

host similar winds from their inner regions, escape pre- 
vents electrons from fully cooling, and our strong wind 
variant would imply that a would decrease to ~ 0.75— 0.8 
at these densities. In our fast diffusive escape model, the 
electrons are similarly prevented from fully cooling, and 
a is slightly reduced in normal galaxies to ~ 0.85 — 0.90. 
However, the efficient escape in these models tends to 
break the FRC. 

Although calorimeter theory often is said to produce 
too high a, we consistently find that a is relatively 
low for starbursts (Figures [2] and ITTj) . The spectral 
slope at 1.4 GHz ranges from 0.7 for weak starbursts 
(S 5 = 0.1 g cm~ 2 ) to 0.5 for extreme starbursts (S g 
1 g cm~ 2 ) in our standard model. In our models, 
the high densities in the starbursts (relative to the low- 
density radio disk of the Milky Way) cause the flat 
spectra. CR electrons and positrons experience severe 
coo ling by bremsstrahlung and ionization, lowering a 
(cf. iThompson et al.l 12006 ). Extre me starbursts are in 
fact observed to hav e fiat spe ctra (ICondon et al.l 119911 : 
IClemens et ail I2008D . though ICondon et al.l (|1991[ ) at- 
tribute the fiat spectra to free-free absorption and ar- 
gues that the intrinsic a is 0.7. We also note that mod- 
els that include the FIR optical depth in U p h predict 
steeper spectra, since IC losses are more effective: our 
kfir = 10 cm 2 g _1 model (Section IA.4j) implies that 
a w 0.65 in starbursts. 

As an example of the power and limitations of our ap- 
proach, we show the predicted synchrotron radio spectra 
of the starbursts in M82, NGC 253, and Arp 220 of our 
fiducial model in Figure Q21 We calculate the radio emis- 
sion using the Schmidt law, and assumi ng a disk geome- 
try wi th the radius of the starburst fromrThompso n~et al.l 
(2006) and scale height h — 100 pc. Our fiducial model 
underpredicts the radio emission of M82 by a factor of 
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Fig. 13. — The fractional contributions of n° decay (solid), 
bremsstrahlung (dotted), and IC (dashed) to the 7-ray flux for 
S 9 = 0.001 g cm' 2 (gray) and £ 9 = 10 g cm -2 (black). The drop 
in bremsstrahlung and IC past 100 GeV comes from the 500 GeV 
cutoff in our electron spectra. 

~ 2 and overpredicts the radio emission of Arp 220 by a 
factor of ~ 4. This is caused by scatter in the Schmidt 
law and the FRC, which our models do not currently ac- 
count for. However, if we normalize the radio spectra to 
the observed 1.4 GHz (dashed line in Figure ITS]) , we find 
that our models predict the radio spectra surprisingly 
well. The spectra of NGC 253 and Arp 220 are slightly 
flatter than predicted, which is probably due to free-free 
absorption. A variant with winds and B oc p a similarly 
predicts the spectral shape, although not the normaliza- 
tion (gray lines in Figure [12]). While the fiducial model is 
no replacement for individual models of galaxies, which 
predict the correct normalization of the radio spectra and 
model the thermal emission and absorption, Figure [12] 
demonstrates that the GHz radio spectra of starbursts 
in general can be understood well in terms of the high- 
Ti g conspiracy. 

5.5. The j-Ray (and Neutrino) Luminosities of 
Starbursts 

Starburst galaxies are predicted to be strong sources 
of 7-rays, observable with Fermi and Very High En- 
ergy (VHE) tele scopes. Previous studies have consid - 
ered NGC 253 (Ipomingo-Santamaria fc Torres! [2005). 
M82 (iPersicet al.l l2008Ude Cea del Pozo et al.l [20091 ). 
Arp 220 (lTorresll2004D . and the diffuse 7- ray background 
(|Thompson. Quataert. &; Waxmanll2007| ). Several star- 
burst galaxies have already been observed in VHE 7- 
rays to search for the emission. Until recently, only up- 
per limits were available on their 7-ray emission (e.g., 
lAharonian et al.ll2005t lAlbert et al.ll2007l) . However, de- 
tections of NGC 253 an d M82 have now been announced 
with VHE telesc opes (lAcciari et al.l l2009t lAcero et all 



120091 ) and Fermi (|Abdo et al.ll2010H . 

Pionic 7-rays come from CR protons in the ISM of the 
starbursts. Since our explanation of the FRC requires 
that secondary electrons and positrons contribute to the 
radio emission, the 7-ray luminosities of starbursts are a 
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useful test of the high-E g conspiracy. 

We calculate the 7-ray flux 18 from secondary 7r° de- 
cay for M31, NGC 253, M82, and Arp 220 in Tabled 
as a check on o ur models. We use the Schmidt law 
and th e E g from iKennicuttl ()1998f ) and iThompson et al.1 
(2006) to calculate the emissivities of gamma rays for 
these systems, whi ch we then multiply by volume (from 



the radii given in IThompson et al 



piy oy 
1 120061 



and the scale 



heights in Section I2.2.2[) to get total luminosities to be 
converted to fluxes. Since we are using approximate re- 
lations such as the Schmidt law, our models will be less 
accurate than more detailed models of individudal galax- 
ies, and the predicted 7-ray luminosities are rough es- 
timates only. These models are not meant to replace 
individual models of starburst galaxies. The main ad- 
vantage of our approach is only that we consider star- 
bursts like M82 and NGC 253 in the broad context of all 
star-forming galaxies spanning the range between normal 
galaxies and ULIRGs; our models are necessarily more 
qualitative than more specific predictions. 

Inelastic proton-proton collisions will also create neu- 
trinos and antineutrinos. The total neutrino (y + v) flux 
is approximately equal to the ir ° 7-ray flux at energies 
E > m^c 2 « 140 MeV (|Steckerlll979l lLoeb fc Waxmanl 
[2006T) . Although we do not calculate the neutrino flux 
directly, we note that the values listed in Table [2] would 
also be good estimates for the neutrino fluxes, summed 
over all flavors and including both neutrinos and antineu- 
trinos. 

Bremsstrahlung and IC emission also are expected to 
contribute to the gamma-ray luminosities, especially at 
low energies. We calculate the bremsstrahlung spectrum 
for M31, NGC 253, M82, and Arp 220 with our stan- 
dard parameters. Both our standard model and our fidu- 
cial wind model imply that in starbursts bremsstrahlung 
emission equals the total pion emission at 100 MeV and 
decreases at higher energies (see Figure[l3]). In less dense 
galaxies, bremsstrahlung grows in importance, but is still 
a minority contributor above 100 MeV. The high energy 
fall-off for bremsstrahlung comes from the steepness of 
the electron and positron spectra relative to the proton 
spectra. About half of the energy in the bremsstrahlung 
emission is below 100 MeV, because the electron spec- 
trum steepens above 100 MeV. 

The IC emission, when integrated over energy, is less 
than the bremsstrahlung or pion 7-ray emission (Fig- 
ure [4]) . An IC gamma-ray spectrum would require an in- 
cident spectrum including CMB, dust, and stellar emis- 
sion. To get a feel for the IC emission, we model the 
background emission as three blackbodies: the CMB, a 
dust component (20 K in normal galaxies and 50 K in 
starbursts), and a direct stellar component (10000 K). 
The dust component and the stellar component have a 
total energy density of U p h : * (eq. [9]) and are scaled ac- 
cording the UV optical depth tqv (see Section |2~3| . 

We find that IC is smaller than bremsstrahlung for 
energies above about 1 MeV, and smaller than pion 7- 
rays above about 50-80 MeV. At low energies, dust emis- 
sion dominates the IC emission in galaxies with E g > 

18 We do not include any o ptical depth to 7-rays in our calcula- 
tions. However, Torres (2004) found that Arp 220 was opaque to 
7-rays only at energies above 1 TeV, and this should also be true 
for galaxies with a lower surface density. 
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Fig. 14. — The fractional contributions of the CMB, dust emis- 
sion, and starlight to the IC emission of galaxies. For starbursts 
with £ 9 > 0.1 g cm -2 , upscattered emission from dust dominates 
at all shown energies. The drop in the contribution of the CMB 
near ~ 1 GeV and dust past ~ 10 GeV is an artifact of our 500 
GeV cutoff in the electron and positron spectra. 

0.01 g cm -2 (see Figure IT4"f . Upscattered UV emission 
from stars only dominates at high energies, and the CMB 
dominates the low energy IC emission in low surface den- 
sity galaxies. Our predicted spectra for starbursts have a 
precipitous fall-off in IC emission for the starbursts past 
10 GeV, because there are no electrons above our 500 
GeV cutoff to boost FIR photons to higher energies. In 
low surface density normal galaxies, the CMB continues 
to provide the photons for most of the IC emission, so 
the drop is at ~ 1 GeV, with UV starlight providing 
higher energy photons (see also Figure [7]) . Most of the 
energy in the IC emission is at low energies, with more 
than half of the upscattered IC photons having less than 
5 MeV. This is because the electron spectrum steepens 
for E > 100 MeV (7 « 200) and incident photons are 
upscattered in energy by a factor of about 7 2 ; an inci- 
dent 10 eV photon would typically only be boosted to a 
few hundred keV. 

Considering the uncertainties and approximations in 
our approach, our models are in loose order-of-magnitude 
agreement with more sophisticated models (Table 
The high-energy pionic 7-ray spectra (> 1 GeV) are 
largely the same as previous models for M82 and 
Arp 220, although our M82 models are near the low 
range of the predictions of lde Cea del Pozo et al.l (|2009f) . 
We predict low er fluxes for NGC 253 than previ- 
ous models by IDomingo-Santamaria fc Torres! (2005) 
and iRephaeli et al.l (|2010j) . by a factor of ~ 4 — 13. 
In particular, we predict an integrated flux of F 7 (> 
100 MeV) « 4 x 10~ 9 for NGC 253. This is sub- 
stantially smaller than the prediction s of 2.3 x 10~ 8 by 
IDomingo-Santamaria fe Torres! (|2005! ) and 1.8±J | x 10" 8 
bv IRephaeli et al.l (|2010ft . 

Interestingly, our predictions for M82 and NGC 253 
are comparable to the Fermi and VHE detection s 
(|Acciari et al.ll2009t lAcero et al.ll2009t lAbdo et alj2010D . 
However, this agreement is caused by a fortuitous cancel- 
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lation of factors. In Table [2] we show that the predicted 
SN rates for M82 and NGC 253 using the Schmidt law 
are very smal l . If we combine the IR luminosities of 
ISanders et"all ()2003[ ) with equations [2] and H we find SN 
rates of 0.065 yr" 1 for M82 and 0.039 yr" 1 for NGC 253. 
Only half of the F IR luminosity of N GC 253 comes from 
its starburst core (|Melo et al.ll2002D . so its starburst has 
a supernova rate of ~ 0.019 yr -1 . By contrast, our model 
using the Schmidt law implies SN rates that are smaller 
by a factor of - 3 for M82 and - 2 for NGC 253's star- 
burst core. 19 If we scale the 7-ray fluxes to these SN 
rates, then the 7-ray fluxes are near or somewhat above 
the upper ranges of previous models. Furthermore, our 
rescaled fluxes for M82 and NGC 253 are then about 
twice as high as observed. 

Other differences with the models arise because 
we also use different distances to the starbursts (we 
use 3.5 Mpc for NGC 253 instea d o f 2.5 Mpc, as 
Domingo- Santamaria fc Torresll2005l and iRephaeli et al.l 
20101 did), but the other models fit the observed radio 
emission so a greater distance would be fit with a greater 
luminosity in these models. We also use different low- 
energy energy spectra (we simply use E~ p instead of K~ p 
or q~ p where q is momentum), which will tend to under- 
estimate the low energy CR proton spectrum, although 
the higher energy CR proton spectrum will be largely 
the same. We again emphasize that our current generic 
models cannot replace existing models, but are meant 
as a demonstration of principle for the broad range of 
star-forming galaxies. 

Our predicted fluxes do not change significantly if 
we consider models with winds and B cx p a (see Sec- 
tion though the TeV fluxes of starbursts are higher 
with this variant because we use p — 2.2 instead of 2.3. 
The 7-ray fluxes also provide good tests of the high-E g 
conspiracy in our models. The ir° 7-ray fluxes we predict 
for starbursts are mainly determined by proton calorime- 
try, the fraction of electron power lost to synchrotron, 
and the Milky Way CR proton normalization. Note that 
our fluxes are several times greater th a n thos e predicted 
by iThompson. Quataert. fc Waxmar] {2007), who as- 
sumed proton calorimetry but did not take into account 
non-synchrotron losses. The significant bremsstrahlung, 
ionization, and IC losses in our model requires more (sec- 
ondary) electrons and positrons to get the same radio 
emission, in turn requiring more protons. 

Data from Fermi and VHE telescopes can distinguish 
these scenarios. Proton calorimetry implies a hard E~ 22 
7-ray spectrum instead of the Galactic E~ 2J spectrum: 
proton calorimetry increases the high-energy "/-ray emis- 
sion. The detections of starburst galaxies with VHE 
telescopes support a hard spectrum and proton calorime- 
try. Note that M31 is much fainter than the starbursts 
in VHE 7-rays, because of its steeper escape-dominated 
CR proton spectrum, though the flux of E < 1 GeV 7- 
rays from M31 is similar to that from nearby starbursts. 
Bremsstrahlung and IC 7-rays overwhelm ir° 7-rays at 
E < 100 MeV in our model, as seen in Figure [T"3l non- 
synchrotron losses increase the low- energy j-ray emis- 

19 If we scale to the 1.4 GHz radio luminosity from 
IWilliam s & Bower (2010) instead of the TIR luminosity, we find 
that M82 must be scaled up by ~ 1.8 and NGC 253 by ~ 1.3 in 
our fiducial model. The agreement is better still in our fiducial 
model with winds and B oc p a . 
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Fig. 15. — The importance of magnetic, radiation, and CR pres- 
sures compared to the hydrostatic pressure needed to support a 
galactic disk. The hydrostatic pressure needed to support the gas 
alone is TrGSg. In low-density galaxies, the mass of the stars im- 
plies that Phydro = 107rG£g (see the discussion in Section 15. 6|l , 
The cosmic ray energy density does not increase as quickly as ra- 
diation and magnetic field energy densities in starburst galaxies. 
None of the three components provides enough pressure to support 
starburst galaxies. 

sion. Fermi detection of this low energy emission would 
support the importance of non-synchrotron cooling. Fi- 
nally, a different proton normalization simply changes 
the amount of both CR protons and secondary electrons 
and positrons: high proton normalization increases j-ray 
emission at all energies, without changing the spectrum. 

The current Fermi and VHE detections of M82 and 
NGC 253 are somewhat ambiguous, because these star- 
bursts are relatively wea k and there is no spectral infor- 
mation at 100 MeV yet ()Abdo et al.ll2010D . The implied 
GeV-to-TeV spectral slopes are ~ 2.2 — 2.3, which is 
consistent with proton calorimetry (|Acciari 
lAcero et aTll2009l; lAbdo et al.ll2010n . However, the fluxes 
are lower than our predicted fluxes scaled to the IR lu- 
minosities of these galaxies. This can imply that either 
proton calorimetry is weaker, or the high-E g conspiracy 
is weaker, particularly in NGC 253. We note that sev- 
eral groups have estimated E g ~ 0.1 — 0.2 g cm~ 2 for 
these starbursts, so that the observed winds could be 
sufficient to break proton calorimetry. More data and 
more sophisticated modeling are needed to fully under- 
stand the implications of these 7-ray observations. Fu- 
ture detections of ULIRGs, which are more likely to be 
proton calorimeters, would be particularly helpful in un- 
derstanding whether there is a high-S g conspiracy. 

5.6. The Dynamical Importance of Cosmic Ray 
Pressure 

Near the Solar Circle in the Milky Way, the CR en- 
ergy density approximately equals magnetic field en- 
ergy density, gas pressure, and radiation energy den- 
sity. Their pressure is also comparable to the pres- 
sure needed to support the Milky Way hydrostatically, 
-Phydro ~ TrGEgEtot ~ IOttGY, 2 ,, where S to t ~ 10S 9 is the 
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surface density of all matter i n the Galactic disk. Extrap- 
olating from t he Milk y Way (IChevalier fc Franssonlll984t 
lEverett et aHl2008ft . ISocrates et al.l (|2008l ) hypothesize 
that CRs continue to provide significan t pressure sup- 
port a nd that they drive strong winds. Jubclgas et al. 
(2008) have also explored the dynamical importance of 
CRs in galaxies, and conclude that they are not impor- 
tant for starbursts. 

We show in Figure [15] the pressure from magnetic 
fields, radiation, and CRs in our standard model, com- 
pared to the hydrostatic pressure needed to support a 
galactic disk. Magnetic fields and radiation (including 
FIR light ) remain compara ble as S ff increases, as pre- 
dicted bv iThompson et al.l (|2006l) . For the inner Milky 
Way, the predicted CR pressure is 2.7 x 10~ 12 erg cm~ 3 , 
within a factor of 2 of the derived best-fit CR pres- 
sure m lEverett et ail ([20081 ). CR pressure remains m 
rough equipartition with magnetic field and radiation 
pressure until the weak starbursts, but then increases 
much more slowly. As CR pressure is mainly provided 
by protons, the failure of CR pressure is caused by pion 
losses: starbursts are proton calorimeters, converting 
most of the CR proton energy into gamma rays and neu- 
trinos that escape the system. Our results are consistent 
wit h the low aver a ge CR pressure in starbursts found 
by iJubelgas et al.1 (|2008| ) , though the CR contribution 
may increase near t he starburst edges (Appendix C of 
ISocrates et ai1l2008h . 

6. SUMMARY AND FUTURE IMPROVEMENTS 

We model the FRC across the range 0.001 g cm~ 2 < 
Eg < 10 g cm -2 , from normal spirals to the densest 
starbursts. The correlation holds in several scenarios de- 
scribed in Appendix |A"1 We find that: 

• We are able to reproduce a linear FRC (Figure [1]) 
consistent with both local and integrated Galactic 
constraints on the energy in CR protons. We find 
that £ « 0.021 of an SN's energy goes into CR 
electrons and rj 0.1 goes into CR protons when 
p w 2.3 and using an E~ p spectrum. 

• Starburst galaxies (S g > 0.1 g cm~ 2 ) are UV, 
electron, and proton calorimeters for most possible 
scenarios. In our standard model, normal galax- 
ies with E ff = 0.01 g cm -2 are UV and electron 
calorimeters, but they are not proton calorimeters 
with only ~ 5% — 15% of CR proton energy going 
into pion losses. 

• The FRC is caused by calorimetry combined with 
two conspiracies operating in different density 
regimes. At low Eg, decreasing electron calorime- 
try causes lower radio emission, but is balanced 
by decreasing UV opacity, which causes lower IR 
emission (Figure [I). At high E ff , bremsstrahhmg, 
ionization, and IC losses decrease the synchrotron 
radio emission, while the appearance of secondaries 
and the effects of B on vc increase the radio emis- 
sion (Figure [5]). 

• The magnetic field strength scales as B oc Eg 6-0 7 , 
implying B w 1 — 2 mG in extreme starbursts with 
E g = 10 g cm -2 . Magnetic fields are significantly 



below equipartition with respect to gravity in star- 
bursts. 

• The CR pressure remains in equipartition with ra- 
diation and magnetic field pressure for galaxies 
with Eg < 0.1 g cm -2 . In starbursts, the CR pres- 
sure is significantly below equipartition, because of 
pion losses (long dashed line in Figure [T5|) . 

• Despite the short synchrotron and IC cooling 
timescales, our models reproduce the observed flat- 
tened radio spectra of starbursts (Figure [2]), be- 
cause of the strong bremsstrahlung and ionization 
cooling in these galaxies. 

• Our models predict that FRCs exist at frequencies 
other than 1.4 GHz, though with increased non- 
linearity (Figure fT0| . 

• Our predictions for the 7-ray emission from M82 
and NGC 253 are within an order of magnitude of 
the Fermi, VERITAS, and HESS detections. How- 
ever, these models assume the Schmidt law holds 
exactly for these starbursts. If we normalize our 
models' IR emission to the observed IR emission, 
which should scale as star formation and CR injec- 
tion power, we find that our 7-ray predictions are 
~ 1.5 times higher than observations of M82 and 
NGC 253, possibly because of strong winds in these 
starbursts. Our predictions Arp 220 are roughly in 
line with previous theoretical models, considering 
the approximations we make (Table [2]) . Full un- 
derstanding of the 7-ray fluxes of these individual 
galaxies probably requires more refined models. 

Our models still have several unresolved issues. We 
have trouble matching the spectral slope a to observa- 
tions of normal galaxies: we predict spectra that are 
too steep. A possible solution may be the presence of 
a wind lowering the escap e time. Howev e r the addition 
of a wind, as observed bv lEverett et al.1 (|2008f ). to our 
models of the Milky Way tends to break the FRC. It is 
also possible that stronger diffusive escape is present in 
the radio halos of normal galaxies than we used in our 
models, because the normal galaxy radio scale hei ght i s 
typically less than the CR scale height (see Section lA.3p . 

Our one-zone models include CR cooling processes and 
escape through diffusion (winds are considered in Ap- 
pendix and can test a variety of parameterizations 
for the environment CRs travel through. Not every issue 
was considered in this paper, though. A natural ques- 
tion would be how robust the FRC is to scatter in the 
properties of the host galaxy environments. F or example, 
the S chmidt law has a scatter of > 0.3 dex (|Kennicuttl 
1998), compar able to the FRC 's own scatter of about 
0.26 dex (e.g.. lYun et~aT1l2001h . It is also unlikely that 
the magnetic field exactly scales as p a or S°, or that the 
overdensity / of ISM gas that CRs travel through would 
be exactly the same from galaxy to galaxy. 

We focused on star formation and the CRs it pro- 
duces in our models of the FRC. However both radio and 
FIR emission have other sources. Star formation drives 
thermal radio emission, whi ch is important at high fre- 
quencies (e.g., [Condon 1992). Thermal free-free emission 
probably also dominates the radio for very low density 
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galaxies like the Large Magellanic Cloud, where the lu- 
mino sity is low and electrons escape easily ([Hughes et al.l 
2006). In normal galaxies, old stellar populations con- 
tribute a significant amount of FIR light without gener- 
ating CRs. To account for this, we might have to distin- 
guish between a warm component of FIR directly related 
to star formation and a cool FIR component that in- 
cludes old stars, and make predictions fo r the FIR colors 
of galaxies and starbursts (|rIeloulll986l ). 

An obvious modification would be to apply our models 
to higher redshifts. Although the FRC has been mainly 
studied in the low-z universe, there have been several 
recent studies o f high-z star-fo r ming galaxies. Recent 
observations by IVlahakis et al.l (|2007|) have found that 
starbursts become radio brighter at high redshift rel- 
ative to the 5 ~ FRC (see also, e.g.. [Kovacs et all 
l2006HMurphvll200alMTchalowski et al.ll2010D . In general, 
though, calorimeter theory predicts that LTm/^radio 
should not change in starburst galaxies. Other studies 
have found that the FRC holds unchanged at high red- 
shifts (e.g. JAppleton et aLl l2004). At high redshifts, the 
CMB will have a greater energy density, implying greater 
IC losses. However, the CMB will not be important in 
dense galaxies, except at the greatest redshifts (c.f. the 
CMB line and the starlight radiation line in Figure [15)) . 
More important are the morphology changes. Many star- 
bursts at high redshifts are observ ed to be kiloparsecs in 
radius instead of ~ 100 pc (e.g.. [Chapman et al.l [2004t 
iBiggs fc Ivisonl 120081: [Younger et all 120081) . usually with 
moderate £ s but at l east one with a surf ace density com- 
parable to Arp 220 (| Walter et al.ll2009D . In these star- 
bursts, the high-S g conspiracy can be unbalanced, alter- 
in g i^TiR./^radio- We will ex plore these effects in detail 
in lLacki fc Thompson! ([2009D . 

While we assumed that galaxies and starbursts are ho- 
mogeneous, future improvements can be made by using 
simple few-zone models of the ISM. In normal galax- 
ies, the CRs are injected from a gas-rich thin disk, but 
can diffuse within a thicker radio disk containing much 
lower density ionized gas. This is reflected in the ra- 
dio emission in the Milky Way, which has both a thin 
and a t hick disk, the latter providing most of the lumi- 
nosity IMS EUl [BeuemanEitil][i98l). Two -zone 
models can account for these density variations. Even 
within the gas-rich disk, the density can fluctuate wildly 
between the high-density molecular clouds and the low- 
density coronal phase. The magnetic field and density 
also change from spiral arms to interarm regions, as well 
as with distance from the centers of galaxies. Similarly 



in starbursts, most of the gas is believed to be in a phase 
with high de nsity, while most o f the volume is relatively 
low density (|Greve et al.ll2009t ). If the CR populations 
in such phases are not well mixed, bremsstrahlung and 
ionization losses would be weak in a low-density phase, 
but strong in a high-density phase. The low-density 
phase, with more volume, would contain most of the CRs 
but might have a steep spectrum, while the high density 
phase, with more mass, would have fewer CRs but with 
a harder spectrum. However, including these different 
phases would add additional parameters to the models. 
The structure of the ISM phases would have to remain 
generic, because detailed information is only available for 
the Milky Way, but the FRC spans a vast range in star 
formation rate and gas surface density. These parame- 
ters would complicate the conspiracy even further, since 
there would be more parameters to tune. 

Ultimately, abandoning the one-zone (or few-zone) ap- 
proach would be necessary for a full understanding of the 
FRC. Our approach only considers the global properties, 
but the FRC holds locally in galaxies to sub-kiloparsec 
scales. We could address the local properties of the FRC 
by making full diffusion models, similar to the GAL- 
PROP models for the Milky Way, for galaxies across the 
entire range of the FRC. A compl ete theory might have 
to include time evolution as well: iMurphv et al.l (2008) 
found that synchrotron emission is better correlated spa- 
tially with star formation in regions of high star forma- 
tion, possibly because the CR electrons have not yet had 
time to diffuse. Spatial diffusion and time dependence 
would make modeling vastly more complicated, but in- 
cluding them may eventually be worthwhile with future 
improvements in radio and infrared observations. 
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APPENDIX 
VARIANTS 

Our standard model (Sections 121 and l4Tj) does not include a variety of effects that may alter the CR populations in 
star-forming galaxies. Our one-zone models of CR injection, cooling, and escape allow us to efficiently survey many 
scenarios. These include variations in essentially all of the parameterizations of Section [2] We search for models that 
are successful under the observational constraints described in Section 12.31 We test several combinations of these 
effects with sparser grids of models, spanning p — 2.0, 2.2, 2.4 and 2.6, and choose / so that fh = 1.0 kpc and 2.0 kpc. 

We summarize our results in Table [3] Table [3] shows the values of p, f, a, and S that satisfy local-based constraints, 
the integrated Milky Way pion gamma-ray luminosity, and both together (Section I2.3[) . Despite the large number 
of scenarios tried, we find similar successful parameters in those scenarios that worked at all. The FRC in some of 
these variants is shown in Figure [TBI Overall, we did not strongly constrain p or /, with the range in allowed / being 
determined mainly by the observed beryllium isotope ratios at Earth. Specific variants occasionally imposed stricter 
constraints on p. However, our models did place strong limits on the magnetic field energy density (in the form of a; 
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Fig. 16. — The FIR-radio correlation, as reproduced in our standard model and several variants. We plot the values when £ = 0.008 in 
each case. Key: 1 (dash, black) - B oc p a with 300 km s — 1 wind in starbursts (p = 2.2, / = 1.0, 6 = 67, a = 0.5); 2 (long dash-dot, black) 
- B oc p a and strong winds of 175 km s — 1 in E 3 = 0.01 g cm -2 , 600 km s — 1 for starbursts (p = 2.2, / = 1.5, 6 = 45, a = 0.5); 3 (dotted, 
black) - standard model with kfir = 1 cm 2 g — 1 (p = 2.2, / = 1.0, 5 = 45, a = 0.7); 4 (dotted, light gray) — model with kfir = 10 cm 2 g — 1 
(p = 2.2, / = 2.0, <5 = 22, a = 0.8); 5 (solid, black) - standard model (p = 2.3, / = 1.5, <5 = 48, a = 0.7); 6 (long dash short dash, black) - 
fast diffusive escape (p = 2.2, / = 2.0, a = 0.6, S = 45); 7 (long dash, gray) - constant D z , B oc p a , and winds of 300 km s — 1 in starbursts 
(p = 2.2, / = 1.5, <5 = 34, a = 0.6); 8 (dash-dot, gray) - same as (7) but with FIR opacity of ftpiR = 1 cm 2 g -1 . 

eq. fT2|) and the CR energy density (in the form of 5; eq. fT7|) . We found a to be 0.5 - 0.8, depending on whether B was 
parametrized to vary with density or surface density. Table [3] shows the allowed 5, because S decreases by a factor of 
100 as p increases from 2.0 to 2.6 from the normalization issue discussed in Section [57X1 In most cases, S is within the 
range 35 - 100, as expected from the local proton-to-electron ratio p/e of about 100. 

B cx p a 

While our standard model assumes that BaS° (eq. fT2|). the magnetic field in galaxies may vary as p a (jGroves et al.l 
2003). A B tx p 5 scaling is also observed in Galactic molecular clouds (|Crutcherlll999( ). The difference matters most 
when comparing the FRC between starburst and normal galaxies, where the scale height changes by an order of 
magnitude in our models. At the transition between the two regimes, since h decreases by a factor of 10 for the 
starbursts, B in this parametrization jumps by a factor of 10°, strengthening synchrotron cooling (Section 13. 2p . 

This variation on our standard model typically breaks the FRC at the transition between normal galaxies and 
starbursts. Starbursts generally have a iTm/iradio that is several times smaller than normal galaxies, because of the 
dramatic increase in the magnetic field and corresponding synchrotron emission. No model retains the correlation and 
also fulfills the local-based constraints on e + /(e~ + e + ) and p/e, and only one model predicts the integrated Milky 
Way gamma-ray luminosity from pions. Models which do satisfy local constraints on the proton normalization have 
£TiR./£radio that vary by 2.4 at p = 2.6 and more at lower p (> 2.7 at p = 2.2). However, other variants in combination 
with B cx p a turn out to restore the FRC (Section [A. 2[) . 

Winds 

We test the effects of including a wind of 300 km s^ 1 (a relatively "weak" wind) in all starbursts with E 9 > 
0.05 g cm~ 2 , as discussed in Section ^. 2. II (eq. [61). We find that the FRC is broken for all values of the parameters. In 
the best cases, £TiR,/£radio varies by ~ 2.3 over the range in E 9 , while our upper limit to the allowed variation was 
2.0 (eq. fl4|). and this does not consider additional constraints on the proton normalization. The difference comes from 
the E 9 = 0.1 g cm -2 models, because the radio emission is halved as the 1.4 GHz CR electrons escape before cooling: 
essentially, electron calorimetry fails for weak starbursts. While the FRC is broken in models with winds and B cx E g 
because synchrotron cooling is not strong enough, recall that the problem with the B cx p a models (Section IA. 1|) was 
that synchrotron cooling was too strong in low-E g starbursts. This suggests that models with both winds and B cx p a 
might work. 

We test this conjecture by testing cases with winds and B cx p a . We then do find some models that satisfy our 
local constraints, including preserving the FRC. Overall, the derived constraints on p, /, and 6 are similar to our 
standard model (Section 14. lj) . There is one noticeable difference in the allowed parameter space compared to our 
standard models: the models that do work now have a — 0.5 — 0.6 instead of 0.6 — 0.8. However, this still implies 
strong magnetic fields in starbursts, because B does not just depend on E g but also h, where our adopted h is 10 
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times smaller in starbursts than in the Milky Way. If B oc p ' 5 , then the magnetic field strengths are 0.12 mG in 
weak starbursts (£ g = 0.1 g cm -2 ) and 1.2 mG in extreme starbursts (£ g = 10 g cm -2 ). These magnetic fields are 
comparable in strength to those that would be present if B oc £ g ' 7 - A model with p — 2.2, 5 — 67, a = 0.5, / = 1.0, and 
£ = 0.0119 satisfies both local-based constraints and the integrated Milky Way 7-ray luminosity. It also reproduces 
both the CR electron and proton fluxes at Earth above 10 GeV to within a factor of ~ 2. The spectral slope a in 
starbursts is comparable to our standard models, being slightly lower because the CRs escape quicker. 

In this model the high-£ g conspiracy operates in starbursts largely as it does in our fiducial model (see Figure [TT]) . 
However, its onset is more gradual because of the winds in the S 9 = 0.1 g cm~ 2 case, reducing secondary production 
in the weakest starbursts. Furthermore, the low a strengthens the magnetic field in the weakest starbursts, so that 
synchrotron can compete more effectively with bremsstrahlung and secondaries. In dense starbursts, magnetic fields 
are relatively weak and a strong high-E g conspiracy sets the radio luminosity. 

We also test str onger winds in comb ination with B oc p a , motivated by the inference of a wi nd in the inner region s 
of the Milky Way (jEverett et al.ll2008l ). as well as the high wind speeds observed in starbursts (Heckman et al. 2000). 
The new wind speeds are v — 175 km s _1 in the £ g = 0.01 g cm -2 model (comparable to the lEverett et al.ll2008l wind) . 
and v — 600 km s _1 in the starbursts. This variant tends to break the FRC, mainly because CRs escaped too quickly 
in the £ g = 0.01 g cm~ 2 Milky Way analog model. With a wind advecting away CRs at lower surface densities where 
cooling is weaker, most of these models are not sufficiently good electron calorimeters to preserve the FRC. Further 
variants with strong winds are discussed in Section IA.8I 

Other Disk Scale Heights 

Small normal galaxy disk scale heights. In our standard scenario, the scale height jumps down from 1 kpc at 
£ g = 0.01 g cm~ 2 to 100 pc at £ g = 0.1 g cm -2 . The jump can cause discontinuities in the FRC if the parameters are 
not fine-tuned. However, the gas disk of the Milky Way is only 100 pc thick, and there is a thin radio disk associated 
with it. Thus, we try models where h = 100 pc for all £ g . In this run, we consider / values of 0.1 and 0.2 to continue 
to match the inferred ISM density that CRs propagate in the Milky Way from beryllium isotopes. Since h does not 
vary, £ g and p are directly proportional. 

We find that models with p = 2.2 - 2.6, a = 0.5 - 0.6, / = 0.1 - 0.2, and 5 = 34 - 100 satisfy both the integrated 
Galactic pion luminosity, and the local e + /(e + + e~) and p/e constraints. The relatively low a is preferred since / is 
low implying weak bremsstrahlung and ionization. Since there is no jump in bremsstrahlung and ionization losses for 
starbursts because they have the same scale height as normal galaxies, a strong magnetic field would make starbursts 
too radio bright. However, p = 2.4 — 2.6 produces CR spectra that are steeper than observed. Higher p are somewhat 
preferred: there is a window in parameter space with a = 0.5, / = 0.1, and d — 50 — 100 where iTm/^radio varies over 
the range in £ g by 2.26 to 2.75 for p = 2.0 but 1.94 to 2.21 for p = 2.4, so that higher p passes marginally. The reason 
for the relatively high variation in Ltir/ -Iradio is that / is low: starbursts are still proton calorimeters, producing 
secondaries that contribute to the radio; bremsstrahlung and ionization losses, which compete for the energy available 
for radio emission, are much weaker than in the standard model (Section I5.2[) . At high surface densities, secondaries 
become important and can make starbursts too radio bright; they are diluted more at higher p (Section 13. ip , so that 
-^TiR/iradio varies slightly less over the range in S 9 . The high-E g conspiracy is present in starbursts for the models, 
but its onset is more gradual, with non-synchrotron losses growing from a factor of ~ 5 at S g = 0.1 g cm~ 2 to ~ 10 at 
E g = 10 g cm~ 2 (see Figure IT?"!) . The more gradual onset arises because we use low / to match the beryllium isotope 
constraints at Earth, weakening bremsstrahlung and ionization and reducing secondaries in weak starbursts; and low 
a, which increases relative synchrotron strength in weak starbursts and suppresses synchrotron in dense starbursts. 

Furthermore, we find that the CR proton fluxes predicted at Earth at 1, 10, and 100 GeV are about 4-20 times 
higher than observed. The predicted electron fluxes are also about 10 times higher than observed. This can easily be 
understood: if h is shrunk 10 times, the same number of CRs are being injected into a smaller volume, so that their 
number density is higher. However, it is possible that the Earth resides in an atypical region of the Galaxy, in which 
case our integrated constraints alone allow this variant. 

We also c onsider a less extreme v ersion of this model, where h norm — 300 pc, about the height of the Milky Way's thin 
radio disk (jBeuermann et aflll985l ). We use / of 0.3-0.6 to match the beryllium isotope measurements. This variant 
is less restrictive, allowing p > 2.2 for all considered constraints. In all of the allowed models, a — 0.6 where BccEJ, 
midway between the a = 0.5 case preferred when h nmra = 100 pc and the a = 0.7 preferred when h norm = 1 kpc. The 
CR flux predicted at Earth is still about 2-5 times higher than observed at Earth, for both protons and electrons in 
these models. 

Large normal galaxy disk scale heights. While the gas disks of normal galaxies are thin, the scale heights of the CRs 
themselves are estimated to be several kpc (see the discussion in Section 12.2.21) . We therefore considered models with 
h = 2 kpc and h = 4 kpc. When h — 2 kpc, we find that the allowed parameters are similar to our standard value. 
Higher a and lower fh is slightly preferred, because the jump in density between normal galaxies and starbursts is even 
greater than when h noim = 1 kpc; therefore, either higher magnetic fields or lower gas densities are needed to prevent 
bremsstrahlung and ionization from overwhelming synchrotron losses. CR proton fluxes at Earth are somewhat small 
by a factor of ~ 1.25 — 2 when p — 2.2. In these models, the high-E g conspiracy is present at an even greater magnitude 
than in our standard model, with non-synchrotron losses suppressing synchrotron emission by a factor of ~ 15 instead 
of- 10 (see Figure [T7|. 

When h = 4 kpc, no models preserve the FRC. Essentially, since / is constant for all S g in our models, and since 
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it must be large to match the beryllium isotope constraints in the Milky Way, bremsstrahlung and ionization are 
inevitably extremely strong in starbursts. While the density increases drastically from normal galaxies to starbursts 
because of the decreasing scale height, magnetic fields do not suddenly increase if they go as E°, and therefore 
synchrotron cannot properly balance bremsstrahlung and ionization for a linear FRC. 

Optically Thick Galaxies 

Normally, our models assume that the CRs propagated in ISM that was optically thin to FIR light. Then, as stated 
in Section \2. 21 E/ph,* = F*/c. However, in a scattering atmosphere, the photon energy density may actually be greater if 
the environment is embedded in an optically thick region. In that case, t/ p h,* = (1 + r-pm)F/c, where tfir = kfirS s /2 
acts as a midplane scattering optical depth. 

In models with kfir = 1 cm 2 g~\ we find that we are still able to recreate the FRC and match both local 
and integrated Galactic constraints. The parameter space allowed by this scenario is similar to our standard model 
(Section 14. lj) . The increased IC scattering in extreme starbursts rules out a = 0.6, so that the magnetic field energy 
density remains comparable to the increased photon energy density. We also typically recover the CR flux at Earth 
to within a factor of 2 for these models when p = 2.2 or 2.4. 

When kfir = 10 cm 2 g _1 , the FRC does not survive in any of our models. The minimum variation in £Tm/£radio 
is 2.03 when a — 0.8, / = 2.0, 5 < 25. Synchrotron losses need to keep up with IC losses in extreme starbursts, which 
would favor high a. However, a — 0.8 — 0.9 often caused too severe synchrotron losses compared to bremsstrahlung 
and ionization in dense starbursts, and too weak synchrotron losses in low surface density galaxies. Those models that 
nearly preserve an FRC have low 5 and high /, which reduces the number of secondaries and increased bremsstrahlung 
and ionization cooling to compensate for the increased magnetic field strength and keep -^Tm/'-E'radio sufficiently high. 

U B = U ph 

Radiation pressure may drive turbulence in the ISM, until the energy densities in radiation and kinetic motions are 
comparable. The turbulence can, in turn, generate magnetic fields. As a result, it is possible that U p h ~ t/turb ~ Ub 
(|Thompsonll2008h . 

We test models where Ub was forced to equal t/ p h, where the radiation energy density includes both the CMB and 
starlight. In optically thin models (tfir = 0), the fiducial values for p, f, and 5 (given in Section l4~Tj) recreate the FRC 
and match both local and integrated proton constraints (see Table[3]). Some models, generally those with p = 2.2 — 2.4, 
also correctly predict to within a factor of 2 the CR proton flux at Earth at energies E = 1, 10, and 100 GeV, as well 
as the CR electron flux at Earth at E = 10 GeV. 

A few optically thick models with kfir = 1 cm 2 g^ 1 where Ub is forced to C7 p h satisfy the integrated Galactic w° 
luminosity, though only one satisfies local constraints. In these models, the effectiveness of synchrotron increased in high 
surface density galaxies, because Ub equaled the quickly increasing U p h- The FRC generally fails because of a tension 
between intermediate densities and high densities. Bremsstrahlung and ionization are stronger at intermediate densities 
(0.01 g cm~ 2 < Yi g < 1 g cm -2 ) than in the standard model (since Ub is lower than our standard prescription), requiring 
a high secondary fraction to compensate. But when S g = 10 g cm -2 , synchrotron cooling overwhelms bremsstrahlung 
and ionization, making those galaxies too radio bright. The vq effect (Section l5.2[) becomes particularly strong as tfir 
becomes appreciable, since B is rapidly increasing, causing the radio emission to increase further. The model that 
does satisfy all constraints has low p, increasing the secondary abundance but weakening the vq effect. 

Increasing the FIR opacity to Kfir = 10 cm 2 g _1 in these variants completely breaks the FRC in all attempted 
models. Again, the bremsstrahlung and ionization losses are strong at intermediate densities (0.01 g cm~ 2 < E 9 < 
1 g cm -2 ) but small at high densities, because both the magnetic field energy density and the radiation field increase 
sharply at S g ss 1 g cm -2 . Therefore starbursts would appear too radio bright to maintain a linear FRC. 

Fast Diffusive Escape 

We have used a scale height of h = 1 kpc for normal galaxies in most variants, based on the scale heights of radio 
disks. However, the diffusive escape time we use in Equation ([5]) applies to the entire CR scale height of the Milky 
Way, which is of order 2 — 4 kpc (see the discussion in Section 12.2.21) . The escape time from the radio disk itself may 
be significantly shorter. Variants in which the scale height of normal galaxies is increased (Appendix IA.3|) still may 
not properly model the radio disk, because they use the midplane magnetic field strength, possibly overestimating the 
effectiveness of synchrotron losses. 

We consider the effect of faster diffusive escape of by running models in which the diffusive escape time in galaxies 
is shortened by a factor of 4 from the nominal lifetime in equation [5] We find that the FRC tends is broken in almost 
all models, because escape reduces the radio luminosity of the lowest surface density galaxies. The minimum variation 
in LpiR/Xradio hi an Y model is 1.998, barely under our criterion of 2.000; this models also satisfies the integrated 
Galactic pion luminosity (p = 2.0, / = 2.0, a = 0.6, 5 = 35, £ = 0.018). The models with the most-preserved FRC 
tend to have low a of 0.6, weakening the synchrotron emission in high surface density galaxies. On the other hand, 
the spectral slope of normal galaxies in these models is 0.85 — 0.90 when p = 2.2, which is closer to the observed values 
than in our standard model (see Section [5^4]) . The CR flux at Earth in models where iFm/iradio varies by < 2.2 and 
p = 2.2 ranges from ~ 35% to 140% of the observed values (68% — 142% when the local and integrated Galactic proton 
constraints hold). 
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Varying Escape Times 

Constant D z . So far, we have been simply assuming that escape time by diffusion for CRs is the same in all galaxies 
and starbursts. However, in our models, the scale height of starbursts is 10 times smaller than that of normal disk 
galaxies. Another simple assumption would be that the vertical diffusion constant D z at any given energy is constant 
across star-forming galaxies and starbursts. Then, since tdifF = h 2 /D z , the escape time would be a hundred times 
smaller for starbursts. This could break the FRC at the transition between normal galaxies and starbursts. 

To see whether this variation had any effect, we modify the escape time to t^\s(E) = tdiff,Mw(^)(^-/^-Mw) 2 - We ran 
our grid for two cases that had worked previously: the standard model where B cx E g and with no winds, and the 
case with (weak) winds and B oc p a (see Section IA. 2 [) . In standard models with constant D z , the vastly more efficient 
escape of CRs in weak starbursts (S g = 0.1 g cm~ 2 ) broke the FRC. Models with B cx p a and winds did reproduce 

the FRC and were able to satisfy local and integrated constraints. The allowed values for p, /, a, and S are similar to 
those for our models with B cx p a , winds, and constant idiff (c.f. Section [A.2I and Table 13]) ■ The increased magnetic 
field strength in starbursts (from B cx p a , where p jumps up for starbursts) compensates for the decreased diffusive 
escape time to restore the FRC in these models. The models generally predict low CR flux at Earth when p = 2.2, 
with a proton flux of about 50% - 120% of its observed Earth value. 

The models with constant D z weaken the high-E fl conspiracy somewhat, but it remains present (Figure lTTl) . The very 
strong diffusive losses and the winds mean that starbursts with S 9 < 1 g cm -2 are not proton calorimeters. Since B 
scales with p a , synchrotron losses can compete more effectively with bremsstrahlung and ionization in weak starbursts. 
Bremsstrahlung and ionization mainly balance the vq effect in these galaxies, which only lowers £TiR./£radio by a 
factor of ~ 50% from E 9 = 0.001 g cm -2 normal galaxies. In denser starbursts, bremsstrahlung and ionization grow 
in importance, but are balanced by secondaries. Overall the conspiracy presented in Figure [T7] is weaker than in the 
fiducial model, with non-synchrotron losses suppressing synchrotron losses by only ~ 2 — 8 for starbursts, because a is 
relatively high. 

Note that constant D z from normal galaxies to starbursts is not necessarily correct. T he speed CRs can stream 
out of galaxies is expected to be limited to the Alfven speed, va — B/y/Anp (Section 15.11 iKulsrud fc Pearcei ri969). 
For B cx p ' 5 , this comes out to a diffusive escape time of 3.7/iioo x 10 6 yr, where hioo = h/(100 pc). While this is 
shorter than the Milky Way CR escape time (eq. EJ, it is also 10 times longer than the constant D z escape time (see 
also Figure 0]). These considerations suggest that diffusion is weaker in starbursts than in constant D z models, which 
would help pre serve the FRC. 

D z cx p^ 1 ! 3 . iHelou fc Bicavl (|1993|) proposed that the diffusion constant scales as p~ 1 ^ 3 , which would make escape 
less efficient in high-density galaxies. We try models where the escape time equaled its local value (eq. [5|) for the local 
surface density (E g = 0.0025 g cm~ 2 ) and increasing as p 1 ! 3 to account for this effect. Testing this assumption against 
both our standard assumptions (no winds and B cx E g ) and B cx p a with winds, we find similar results to the constant 
D z case, although a and / are more severely constrained. The faster escape time still breaks the FRC in otherwise- 
standard models. When B cx p a , though, the models restore the FRC and satisfy local and integrated co nstra ints. 
Again, the allowed parameters are similar to the case with B cx p a , winds, and constant idiff (cf. Section [A. 2 1 and 
Table [3]). The CR flux at Earth is correct, within about a factor of 2 of the observed values for the energies we 
considered. 

D z cx p~ x ■ Finally, we try a rapidly scaling diffusion constant, motivated by the parametrization in iMurphv et all 
(2008). As before, the escape time is normalized to its local value (eq. [S]) at the local surface density (E g = 
0.0025 g cm -2 ), but now increasing as p. We once again test it against B cx E g with no winds and Bap" with 
winds. This variant fails in both cases to create the FRC, because escape is too efficient in the lowest density galaxies 
(E 9 = 0.001 g cm -2 ). 

Multiple Effects and Other Variants 

Combinations with weak winds. We finally consider scenarios that combine most of the previous variants. We include 
starburst winds, the FIR optical depth, B cx p a , and constant D z , D z cx p^ 1 ^ 3 , or D z cx p~ x simu ltaneously. Our 
results are essentially the same as the models with constant or varying D z considered in Section IA. 71 with B cx p a and 
winds. For constant D z , the models reproduce the FRC and satisfy both local and integrated constraints. As before, 
the models predict low CR proton flux at Earth. Some models with D z cx p~ x l 3 also are consistent with the FRC, but 
only only fulfills both local and integrated Galactic proton normalization constraints. The increased photon energy 
density suppresses radio in E 9 = 10 g cm -2 starbursts, thus making it hard to maintain the FRC over the entire range 
in E g . As before models with D z cx p~ x fail to reproduce the FRC. 

Combinations with strong winds. The main problem with the strong wind scenario (Section IA.2|) is the rapid escape 
of the CR electrons when E 9 = 0.01 g cm -2 . More synchrotron is emitted before escape if the CRs have to travel a 
larger distance. We therefore consider models with strong winds and h nmm = 2 kpc, with / scaled to 2.0, 3.0, and 
4.0 to match local isotope measurements. We use a constant diffusion rate D z , and try models with and without FIR 
opacity. In both cases, we are able to satisfy both local and integrated constraints. 
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TABLE 1 
List of Symbols Used 



Symbol Section Standard Definition 

Value 



Derived parameters 



rmrm 

[231 

233T 



2~2~3l 



3.2 



3.3 



2.3 
0.023 
0.12 
5 

48 
0.V 
1.5 



Power law index of the injected spectrum of primary cosmic rays 
Fraction of supernova kinetic energy injected into primary CR electrons 
Fraction of supernova kinetic energy injected into primary CR protons 
?;/£, the ratio of energy in injected protons to injected electrons 

Proton-to-electron ratio at relativistic energy at injection; <5 renormalized to remove p dependence 
Power law scaling of galactic magnetic fields with surface or volume density 
Ratio of density through which CRs propagate to average ISM density 



Other input parameters 



s 9 




m 






Average gas column density 


SsFR 




m 






Star formation rate per unit area 


h 


HI2.2.2I 


1 kpc - 100 pc 


Scale height of CR disk 


E 




m 






Total energy of cosmic ray 


*iife(-E) 


HJ2.2.1I 




Escape time for a particle of energy E from the galaxy, includes both advection and diffusion 


Q(E) 


B1I2H 




Energy spectrum of primary CRs injected into the ISM per unit volume 


HE) 




2 






Energy loss rate per particle (positive for energy loss) 


C 




2.1 






Normalization of the injected energy spectrum of primary CRs 


7 




2.1 






Lorentz factor of cosmic ray, E / (mc 2 ) 


7max 




2.1 




10 6 


Maximum Lorentz factor of a cosmic ray at injection 


K 




f2"T 






Kinetic energy of cosmic ray 


e 




2.1 




3.8 x 10- 4 


Radiative efficiency of stellar population 


E51 




2.1 




1 


Mechanical energy per supernova, in units of 10 51 ergs 


1pl7 




2.1 




1 


Conversion rate between the supernova rate per unit mass Tsn and starlight emissivity e p h 



tdiff(E) 

^adv 

(n) 

"eff 

f^ph.CMB 
«FIR 

B 

KTJV 



E 



2.2.1 



2.2.1 



2.2.3 



2.2.3 



T 



Hi 

l2?X4l 



222]|M 



223]|32 



1231 



500 cm 2 g 



Escape time for particle of energy E from the galaxy by diffusion 

Escape time for particle from the galaxy by advection in a wind 

Average number density of hydrogen, S 9 /(2/i) 

Average hydrogen number density the CRs encounter, / (n) 

Energy density in starlight (UV or reprocessed FIR) 

Starlight energy flux 

Energy density of CMB 

Effective ISM opacity to far-infrared light (FIR) 
ISM magnetic field strength 
Effective ISM opacity to ultraviolet (UV) light 
Critical frequency of synchrotron radiation. 



Output 



N{E) 


2 


Final steady-state spectrum of CRs, calculated per unit volume 


e 


2 3 


Emissivity (here, power per volume), for radiation or a CR loss process 




2.3 


Specific emissivity (emissivity per unit frequency), generally for synchrotron radio emission 


£tir 


2.3 


Total infrared emission from young stars 


^radio 


2.3 


Nonthermal synchrotron radio emission 


<?FIR 


2.3 


Rescaled, observed logarithm of ratio iFlR/^radio 


e+/(e++e~) 


2.3 


Fraction of CR positrons in CR electrons and positrons at Earth, usually at 1 GeV 


p/e 


2.3 


Observed ratio of protons to electrons at Earth, usually at 1 GeV 


L ttO 
dI e (E)/dE 


2.3 


Gamma-ray emission from ir° production 


2.3 


Spectrum of CR electrons at Earth from Milky Way sources 


dI p (E)/dE 


2.3 


Spectrum of CR protons at Earth from Milky Way sources 


a 


2.3 


Power law slope of the observed radio flux, dlogi^/dlog^, usually measured between 1.4 and 4.8 GHz 


V(£) 


IT 


Spectral slope of the final steady-state CR spectrum, dlogN (E) / dlogE 


^■sec/ ^ 


3.1 


Fraction of CR electrons that are secondaries 




5.1 


Fraction of CR proton luminosity going into pion losses 
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TABLE 2 

7T° 7-Ray (and -n-* Neutrino 3 -) Fluxes 



Galaxy log 10 £ 9 R D TgN Integrated n° Photon Flux > E 

100 MeV (Total b ) 1 GeV 10 GeV 100 GeV 300 GeV 1 TeV 

(g cm -2 ) (kpc) (Mpc) (yr -1 ) (cm -2 s _1 ) (cm -2 s _1 ) (cm -2 s _1 ) (cm -2 s _1 ) (cm -2 s _1 ) (cm -2 s _1 ) 



Standard 



Milky Way 


-2.0 


4 


0.008 


0.018 


2.7E-4 


(5.1E-4) 


3.4E-5 


8.5E-7 


1.5E-8 


2.3E-9 


2.8E-10 


M31 


-3.0 


20.9 


0.9 


0.019 


3.7E-9 


(1.2E-8) 


3.9E-10 


8.1E-12 


1.4E-13 


2.0E-14 


2.4E-15 


NGC 253 


-0.33 


0.21 


3.5 


0.011 


2.6E-9 


(4.7E-9) 


5.1E-10 


3.2E-11 


1.6E-12 


3.8E-13 


7.5E-14 


M82 


-0.16 


0.23 


3.6 


0.022 


5.1E-9 


(9.1E-9) 


1.0E-9 


6.3E-11 


3.2E-12 


7.8E-12 


1.6E-13 


Arp 220 (east) 


0.78 


0.12 


76.6 


0.12 


6.4E-11 


(1.1E-10) 


1.3E-11 


8.0E-13 


4.3E-14 


1.1E-14 


2.3E-15 


Arp 220 (west) 


0.94 


0.07 


76.6 


0.071 


3.6E-11 


(6.0E-11) 


7.2E-12 


4.6E-13 


2.4E-14 


6.1E-15 


1.3E-15 


Arp 220 (disk) 


1.08 


0.37 


76.6 


3.1 


1.6E-9 


(2.6E-9) 


3.1E-10 


2.0E-11 


1.1E-12 


2.7E-13 


5.8E-14 



B oc p a and winds 



Milky Way 


-2.0 


4 


0.008 


0.018 


2.5E-4 


(4.4E-4) 


3.3E-5 


9.6E-7 


2.1E-8 


3.5E-9 


4.8E-10 


M31 


-3.0 


20.9 


0.9 


0.019 


3.0E-9 


(7.7E-9) 


3.5E-10 


8.8E-12 


1.9E-13 


3.0E-14 


4.2E-15 


NGC 253 


-0.33 


0.21 


3.5 


0.011 


2.4E-9 


(4.1E-9) 


5.2E-10 


4.0E-11 


2.5E-12 


6.5E-13 


1.4E-13 


M82 


-0.16 


0.23 


3.6 


0.022 


5.2E-9 


(8.9E-9) 


1.1E-9 


8.7E-11 


5.5E-12 


1.5E-12 


3.3E-13 


Arp 220 (east) 


0.78 


0.12 


76.6 


0.12 


7.9E-11 


(1.4E-10) 


1.7E-11 


1.4E-12 


9.2E-14 


2.6E-14 


6.1E-15 


Arp 220 (west) 


0.94 


0.07 


76.6 


0.071 


4.6E-11 


(7.9E-11) 


1.0E-11 


7.9E-13 


5.3E-14 


1.5E-14 


3.6E-15 


Arp 220 (disk) 


1.08 


0.37 


76.6 


3.1 


2.0E-9 


(3.5E-9) 


4.4E-10 


3.5E-11 


2.3E-12 


6.6E-13 


1.6E-13 



a Although wc calculate the ir° 7-ray spectrum explicitly, we do not perform a similar calculation for neutrinos. However, the neutrino flux f rom tt* 
decay at energies much highe r than m^-c 2 140 MeV is approximately equal to the 7-ray flux, if antineutrinos and all flavors are included llSteckerl 
IT97alLoeb fc Waxman]|2006ly 

b The flux in parentheses includes bremsstrahlung and IC 7-rays as well as pionic emission. 

c For simplicity, we treat the Milky Way as a point source at the Galactic Center, and consider only its inner regions. 
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TABLE 3 

Successful Models 



§ 


B 


^wind (E^min) 

km s _1 (g cm -2 ) 


KFIR 

(cm 2 g — 1 ) 


ft- norm 
(PC) 


tdiS 


Constraints' 5 


V 


Allowed Values 
/ 


8 


ED 


B oc S° 







1000 


Const, tdiff 


L 
G 
C 


2.0 - 2.6 
2.0 - 2.6 
2.0 — 2.6 


1.0 - 2.0 
1.0 - 2.0 
1.0 - 2.0 


0.6 - 0.7 
0.6 - 0.8 
0.6 - 0.7 


34 - 100 
10 - 152 
34 — 100 


IA.il 


B cc p a 







1000 


Const. t c iiff 


L 
G 






(A 











it, 






V 


^ 


B ocT,g 


300 (0.05) 





1000 


Const, t^iti 


L 

G 
\j 




fA 

















V) 


IA.2I 


B cc p a 


300 (0.05) 





1000 


Const, tdiff 


L 

G 


2.0 - 2.6 
2.0 - 2.6 

•on OA 
z.xj — z.o 


1.0 - 2.0 
1.0 - 2.0 
1.0 — 2.0 


0.5 - 0.6 
0.5 - 0.6 
0.5 — 0.6 


34 - 100 
15 - 100 

tA 1 nn 
o4 — 1UU 


IA.2I 


B cc p a 


175 (0.01) 
600 (0.05) 





1000 


Const, ^iff 


L 

G 

Q 






SI 




171 

V 












B a E° 







100 


Const, tdig 


L 
G 


2.2 - 2.6 
2.0 - 2.6 
2.2 — 2.6 


0.1 - 0.2 
0.1 - 0.2 
0.1 — 0.2 


0.5 - 0.6 
0.5 - 0.6 
0.5 — 0.6 


34 - 100 
15 - 100 
34 — 100 


1A.31 









300 


Const. t<jig 


L 
G 
c 


2.0 - 2.6 

2.0 - 2.6 
on 9 a 


0.3 - 0.6 

0.3 - 0.6 

oq n p. 
U.o — U.O 


0.6 
0.6 - 0.7 

U.O 


34 - 100 
15 - 100 











2000 


Const, tdig 


L 

G 
r 1 


2.0 - 2.4 
2.0 - 2.4 

on OA 


2.0 
2.0 
2.0 


0.7 
0.7 
0.7 


35 - 67 
35 - 67 

00 — Of 


1A.31 


B«S; 







4000 


Const, tdiff 


L 
G 
C 






(A 
V 

















rh 
V 


\M 


BocE« 




1 


1000 


Const. t c iiff 


L 
G 
C 


2.0 - 2.6 
2.0 - 2.6 
2.0 — 2.6 


1.0 - 2.0 
1.0 - 2.0 
1.0 - 2.0 


0.7 
0.7 - 0.8 
0.7 


34 - 91 
10 - 91 
34 — 91 


IA.4I 






10 


1000 


Const, t^iff 


L 

G 
C 

























IA.5I 


y B = u ph 







1000 


Const, tdig 


L 

G 
C 


2.0 - 2.6 
2.0 - 2.6 
2.0 — 2.6 


1.0 - 2.0 
1.0 - 1.5 
1.0 - 1.5 




34 - 100 
34 - 100 
34 — 100 


IA.5I 


U B = U ph 




1 


1000 


Const, tjifl 


L 

G 


2.0 
2.0 

2.0 


1.5 
1.0 - 1.5 
1.5 




50 
20 - 50 
50 


IA.5I 


Ub = U ph 




10 


1000 


Const, tjifl 


L 

G 

c 






















IA.6I 


B ccHg 







1000 


Const, tjiff 
(1/4 nominal) 


L 

G 

c 




2.0 






2.0 






0.6 






35.0 



IA.7I 


B ccHg 







1000 


Const. D z 


L 

G 

c 
























ro 


B CCYlg 







1000 


D z cc p- 1 ' 3 


L 

G 

c 











01 
w 




(1) 







IAJI 


B ccHg 







1000 


D z cc p- 1 


L 

G 








U 










IAJI 


B cc p a 


300 (0.05) 





1000 


Const. Dz 


L 

G 


2.2 - 2.6 
2.0 - 2.6 
2.2 — 2.6 


1.0 - 2.0 
1.0 - 2.0 
1.0 — 2.0 


0.5-0.6 
0.5-0.6 
q 5 _ 0.6 


34 - 100 

10 - 91 
34 — 91 


lA~7l 


B cc p a 


300 (0.05) 





1000 


Dz cc p" 1 / 3 


L 

G 


2.0 - 2.6 
2.0 - 2.4 
2 2 4 


1.0 - 2.0 
in 9 n 

1.0 - 1.5 


0.5 
u.o — u.o 

0.5 


34 - 100 
10 - 67 


|Aj| 


B cc p a 


300 (0.05) 





1000 


D z cc p- 1 


L 

G 
C 























IA.Hl 


B cc p u 


300 (0.05) 


1 


1000 


Const. D z 


L 

G 
C 


2.0 - 2.6 
2.0 - 2.6 
2.0 - 2.6 


1.0 - 2.0 
1.0 - 2.0 
1.0 - 2.0 


0.5 - 0.6 
0.5 - 0.6 
0.5 - 0.6 


34 - 100 
10 - 91 
34 - 91 


IA.8I 


B cc p a 


300 (0.05) 


1 


1000 


D z cc p" 1 / 3 


L 

G 
C 


2.6 
2.0 - 2.6 
2.6 


1.0 - 1.5 
1.0 - 2.0 
1.0 


0.5 
0.5 - 0.6 
0.5 


91 
10 - 91 
91 




B cc p a 


300 (0.05) 


1 


1000 


D z cc p~ L 


L 

G 
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TABLE 3 — Continued 



§ B 


^wind (^gjinin) 


^FIR ^norm 




Constraints' 5 




Allowed Values 








km s _1 (g cm -2 ) 


(cm 2 g" 1 ) (pc) 






P 


/ 


a 




6 










C 
















IA.8I Bap" 


175 (0.01) 


2000 


Const. D z 


L 


2.2 


3.0 - 4.0 


0.6 


34 


- 45 




600 (0.05) 






G 


2.0 - 2.2 


3.0-4.0 


0.6 


20 


- 45 








C 


2.2 


3.0 - 4.0 


0.6 


34 


-45 


IA.8I B oc p a 


175 (0.01) 


1 2000 


Const. D z 


L 


2.2,2.6 


2.0 - 4.0 


0.6,0.5 


34 


- 91 


600 (0.05) 






G 


2.0 - 2.6 


2.0 - 4.0 


0.5 - 0.6 


20 


- 91 










C 


2.2,2.6 


2.0 - 4.0 


0.6,0.5 


34 


- 91 


a Models without a 


wind have a ■ ■ ■ entry. 


Otherwise, there is a 


wind of speed v w 


i n d in models with 


a surface density S g of at least S g , 




In 



the strong wind variants with more than one £ g)m j n listed, the wind speed is that with the greatest S gt min less than S g for each model. 

The FIR-radio correlation must always hold. Additional constraints on the proton normalization: L — local constraints (e~*"/(e~ +e ) and p/e). 
G — integrated constraint (Milky Way gamma-ray luminosity from 7T° production). C — both local and Galactic constraints. 

C Variants with entries could not satisfy the given constraints. Models where the magnetic field has no free parameters have a • ■ ■ entry in 



the a column. 

^ Since the scale height does not vary, D z is constant and B oc p a for these models. 



